WELL POSEDNESS FOR THE MOTION OF A 
COMPRESSIBLE LIQUID WITH FREE SURFACE BOUNDARY. 



Hans Lindblad 

ABSTRACT. We study the motion of a compressible perfect liquid body in vacuum. This can be thought of 
as a model for the motion of the ocean or a star. The free surface moves with the velocity of the liquid and 
the pressure vanishes on the free surface. This leads to a free boundary problem for Euler's equations, where 
the regularity of the boundary enters to highest order. We prove local existence in Sobolev spaces assuming 
a "physical condition", related to the fact that the pressure of a fluid has to be positive. 



1. Introduction 

We consider Euler's equations 



(1.1) p(d t + V k d k )v 3 +djp = 0, j = l,...,n in V, where d i = d/dx i , 
describing the motion of a perfect compressible fluid body in vacuum: 

(1.2) (d t + V k d k )p + pdivV = 0, dwV = d k V k in V, 

where V k = 5 k% Vi = v k and we use the summation convention over repeated upper and lower indices. 



Here the velocity V = (V 1 , V n ), the density p and the domain V = Uo<t<T{t} X V t , V t C W 1 are 
to be determined. The pressure p = p(p) is assumed to be a given strictly increasing smooth function 
of the density. The boundary dT) t moves with the velocity of the fluid particles at the boundary. The 
fluid body moves in vacuum so the pressure vanishes in the exterior and hence on the boundary. We 
therefore also require the boundary conditions on dT> = Uo<t<T{^} x dT>t'- 

(1.3) {d t + V k d k )\ av €T(&D), 

(1.4) p = 0, on dV. 

Constant pressure on the boundary leads to energy conservation and it is needed for the linearized 
equations to be well posed. Since the pressure is assumed to be a strictly increasing function of the 
density we can alternatively think of the density as a function of the pressure and for physical reasons 
this function has to be non negative. Therefore the density has to be a non negative constant ~p on the 
boundary and we will in fact assume that p >0, which is the case of liquid. We hence assume that 

(1.5) p(7>o) = and p\p) > 0, for p > p , where ~p > 

From a physical point of view one can alternatively think of the pressure as a small positive constant 
on the boundary. By thinking of the density as function of the pressure the incompressible case can be 
thought of as the special case of constant density function. 
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The motion of the surface of the ocean is described by the above model. Free boundary problems 
for compressible fluids are also of fundamental importance in astrophysics since they describe stars. 
The model also describes the case of one fluid surrounded by and moving inside another fluid. For large 
massive bodies like stars gravity helps holding it together and for smaller bodies like water drops surface 
tension helps holding it together. Here we neglect the influence of gravity which will just contribute 
with a lower order term and we neglect surface tension which has a regularizing effect. 

Given a bounded domain T> C that is homeomorphic to the unit ball, and initial data Vq and 
Po, we want to find a set V C [0, T] x IR n , a vector field V and a function p, solving (1.1)-(1.4) and 
satisfying the initial conditions 

(1.6) {x; (0,x) eD} = P„, 

(1.7) V = V , p = p on {0}xP . 

In order for the initial-boundary value problem (1.1)-(1.7) to be solvable initial data (1.7) has to satisfy 
certain compatibility conditions at the boundary. By (1.2), (1.4) also implies that divV| ao = 0. We 
must therefore have Po\ gVo = Po an d divVb| a2J = 0- Furthermore, taking the divergence of (1.1) gives 
an equation for (d t + V k dk)divV in terms of only space derivatives of V and p, which leads to further 
compatibility conditions. In general we say that initial data satisfy the compatibility condition of order 
m if there is a formal power series solution in t, of (1.1)-(1.7) (p, V), satisfying 

(1-8) (d t + V k 8 k y (p - P )\ {0}xdT)o =0, j = 0,..,m-l 

Let M be the exterior unit normal to the free surface &D t . Christodoulou[C2] conjectured the initial 
value problem (1.1)-(1.8), is well posed in Sobolev spaces under the assumption 

(1.9) Vatp<-c <0, on dV, where V^ = M i d x i. 

Condition (1.9) is a natural physical condition. It says that the pressure and hence the density is larger 
in the interior than at the boundary. Since we have assumed that the pressure vanishes or is close to 
zero at the boundary this is therefore related to the fact that the pressure of a fluid has to be positive. 

In general it is possible to prove local existence for analytic data for the free interface between two 
fluids. However, this type of problem might be subject to instability in Sobolev norms, in particular 
Rayleigh- Taylor instability, which occurs when a heavier fluid is on top of a lighter fluid. Condition (1.9) 
prevents Rayleigh- Taylor instability from occurring. Indeed, if this condition is violated Rayleigh- Taylor 
instability occurs in a linearized analysis. 

In the irrotational incompressible case the physical condition (1.9) always hold, see [W1,2,CL], and 
[Wl,2] proved local existence in Sobolev spaces in that case. [Wl,2] studied the classical water wave 
problem describing the motion of the surface of the ocean and showed that the water wave is not 
unstable when it turns over. Ebin[El] showed that the general incompressible problem is ill posed in 
Sobolev spaces when the pressure is negative in the interior and the physical condition is not satisfied. 
Ebin[E2] also announced a local existence result for the incompressible problem with surface tension on 
the boundary which has a regularizing effect so (1.9) is not needed then. 

In [CL], together with Christodoulou, we proved a priori bounds in Sobolev spaces in the general 
incompressible case of non vanishing curl, assuming the physical condition (1.9) for the pressure. We 
also showed that the Sobolev norms remain bounded as long as the physical condition hold and the 
second fundamental form of the free surface and the first order derivatives of the velocity are bounded. 
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Usually, existence follows from similar bounds for some iteration scheme, but the bounds in [CL] used 
all the symmetries of the equation and so only hold for modifications that preserve all the symmetries. 
In [LI] we showed existence for the linearized equations and in [L3] we proved local existence for the 
nonlinear incompressible problem with non vanishing curl, assuming that (1.9) holds initially. 

For the corresponding compressible free boundary problem with non-vanishing density on the bound- 
ary, there are however in general no previous existence or well-posedness results. Relativistic versions 
of these problems have been studied in [C1,DN,F,FN,R] but solved only in special cases. The methods 
used for the irrotational incompressible case use that the components of the velocity are harmonic to 
reduce the equations to equations on the boundary and this does not work in the compressible case 
since the divergence is non vanishing and the pressure satisfies a wave equation in the interior. To be 
able to deal with the compressible case one therefore needs to use interior estimates as in [CL,L1]. Let 
us also point out that in nature one expects fluids to be compressible, e.g. water satisfies (1.5), see [CF]. 
For the general relativistic equations there is no special case corresponding to the incompressible case. 
In [L2] we showed existence for the linearized equations in the compressible case and here we prove 
local existence for the nonlinear compressible problem: 

Theorem 1.1. Suppose that p = p{p) is a smooth function satisfying (1.5). Suppose also that initial 
data vq, po and T> are smooth satisfying the compatibility conditions (1.8) to all orders m and V is 
diffeomorphic to the unit ball. Then if the physical condition (1.9) hold at t = there is a T > 
such that (1.1)- (H) and (1.6)- (1.1) has a smooth solution for <t<T. Furthermore (1.9) hold for 
< t < T, with Co replaced by cq/2. 

A few remarks are in order. The existence of smooth solutions implies existence of solutions in 
Sobolev spaces if one has a priori bounds in Sobolev spaces. In the incompressible case we had already 
proven a priori energy bounds in Sobolev spaces spaces in [CL] as well as a continuation result, that 
the solution remains smooth as long as the physical condition is satisfies and the solution is in C 2 . 
Similar bounds in Sobolev spaces hold also in the compressible case. We also remark that there are 
initial data satisfying the compatibility conditions to all orders, see section 16. If only finitely many 
compatibility conditions are satisfied then we get existence in C k for some k. What then is essential is 
that the physical condition hold and this and the existence time only depends on a bound of finitely 
many derivatives of initial data. This makes it possible to construct a sequence of smooth solutions 
converging to a solution in Sobolev norms since we have a uniform lower bound for the existence times, 
in terms of the Sobolev norm. 

A few remarks about the proof are also in order. As in the incompressible case [L3] we will use 
the Nash-Moser technique to prove local existence. However, because of the presence of the boundary 
problem for the wave equation for the enthalpy one has to take as many time derivatives of the equations 
as space derivatives. Therefore one has to use interpolation in time as well and one might just as well 
do smoothing also in time in the application of the Nash-Moser technique although there is no loss 
of regularity in the time direction. Because of this all our constants will depend on a lower bound of 
the time interval. This will make it a bit more delicate since we will also need to choose a small time, 
in order that the physical and coordinate conditions should hold. However, at the same time certain 
estimates are more natural when one includes time derivatives up to the highest order. 

The plan of the paper is as follows. We will assume that the reader is somewhat familiar with the 
notation in [L3]. We will also assume the existence proofs given in [L1,L2,L3] for the inverse of the 
linearized operator so we will only prove improved estimates here. In section 2 we formulate Euler's 
equations in the Lagrangian coordinates and derive the linearized equations in these coordinates. Here 
we also define a modified linearized operator which is easier to first consider. In section 3 we define 
the orthogonal projection onto divergence free vector fields, the normal operator and decompose the 
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equation onto a divergence free part and a wave equation for the divergence. In section 4 we construct 
the families of tangential vector fields, define the modified Lie derivatives with respect to these and 
calculate its commutators with the normal operator and other operators that occur in the linearized 
equation. In section 5 we derive estimates of derivatives of a vector field in terms of the curl the 
divergence and tangential derivatives or the normal operator. In section 6 we give tame estimates for 
the Dirichlet problem, in section 7 we give tame estimates for the wave equation and in section 8 we give 
tame estimates for the divergence free part. Then in section 9 we put these estimates together to get 
tame estimates for the inverse of the modified linearized operator. In section 10 we give estimates of the 
enthalpy in terms of the coordinate. In section 11 we show that the physical and coordinate conditions 
can be satisfied for small times if they hold initially. In section 12 we then get tame estimates for 
the inverse of the linearized operator. In section 13 we give tame estimates for the second variational 
derivative. In section 14 we construct the smoothing operators needed for the Nash-Moser iteration. 
Finally, in section 15 we construct the Nash-Moser iteration that proves local existence of a smooth 
solution. In section 16 we show that one can construct a large class of initial data that satisfy the 
compatibility conditions to all orders. Most of the steps of the proof above works for the pressure any 
smooth strictly increasing functions of the density. However, the estimates for the enthalpy in terms of 
the coordinate simplified if one assume that the pressure is a linear function of the density and we first 
do the proof in this case and then in section 17 give the additional estimates needed for the general 
case. 



2. Lagrangian coordinates and the linearized equation. 



Let us introduce Lagrangian coordinates in which the boundary becomes fixed. Let be a the 
unit ball in R n and let /o : 0, — > Dq be a diffeomorphism. By a theorem in [DM] the volume form 
Ko = det(dfo/dy) can be arbitrarily prescribed up to a multiplicative constant and by a scaling of 
the equations we can also assume that the volume of V is that of the unit ball. Assume that v(t,x), 
p(t,x), (t,x) G V are given satisfying the boundary conditions (1.3)-(1.4). The Lagrangian coordinates 
x = x(t,y) = ft(y) are given by solving 



(2.1) 



dx/dt = V(t,x(t,y)), x(0,y) = f (y), yen 



Then f t : Q — ► V t is a diffeomorphism, and the boundary becomes fixed in the new y coordinates. Let 
us introduce the notation 



(2.2) 



d_ 

di 



y=constant 



d_ 

di 



+ V , 

x=constant dx k 



for the material derivative. The partial derivatives di = djdx 1 can then be expressed in terms of 
partial derivatives d a = d/dy a in the Lagrangian coordinates. We will use letters a,b,c, f to denote 
partial differentiation in the Lagrangian coordinates and i,j,k,... to denote partial differentiation in 
the Eulerian frame. 

In these coordinates Euler's equation (1.1) become 



(2.3) 



pD 2 t Xi + d iP = 0, (t, y) £ [0, T] xfl 



and the continuity equation (1.2) become 



(2.4) D t p + pdivV = 0, (t, y) € [0, T] x Q 
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Here the pressure p = p(p) is assumed to be smooth strictly increasing function of the density p. With 
h, the enthalpy, i.e. h'(p) = p'(p)/p and h = when p = 0, (2.3) becomes 

(2.5) D%Xi + dih = 0, (t, y) € [0, T] xfl 

Since h is a strictly increasing function of p we can solve for p = p(h) as a function of h and with 
e(/i) = In p(h) (2.4) become 

(2.6) D t e{h) + divV = 
Euler's equations are now replaced by 

riii a ri 

(2.7) D t 2 Xi + dih = 0, (t, y) G [0, T] x Jl, where d, 



dx % dy a 



and taking the divergence of (2.7) using that [D t ,di] = —(diV k )dk we obtain a wave equation for the 
enthalpy 



(2.8) D 2 t e(h) - Ah - (0^)0* V* = 0, /i 

Here e(/i) is a given smooth strictly increasing function and 



= 

an 



(2.9) Ah = J2dfh = K- 1 d a { K g ab d b h) where ffab = ^.__ 

and 5 ab is the inverse of the metric g a b and k = det (dx/dy) = \/detg. The initial conditions are 
(2-10) x\ t=Q = f , D t x \t=o = 

(2.11) h\ t=Q = e- 1 Qnp ), D t h\ t=(j = -dwV /e^e" 1 (In Po )) 

where e _1 is the inverse function of e(h). If (x,h) satisfies (2.7)-(2.8) with initial data of the form 

(2.10) -(2.11) then (x,h) also satisfies (2.6). By a theorem in [DM] the volume form can be arbitrarily 
prescribed so we can in fact choose it so kq = det (dfo/dy) = 1/po, in which case e(h) = — lnz-t, since 
this is true when t = and since D t kiK = divF. Hence we are left we two independent initial data, /o 
and vq. 

In order for (2.8) to be solvable we must have the following condition on e(h) and coordinate 
condition: 

(2.12) c^ 1 < e'(h) < ci, 

(2.13) Yl ^ ab \ + \9ab\ < cl \dx/8y\ 2 + \dy/dx\ 2 < c\ 

a.b 

for some constant < c\ < oo. In order for (2.7) to be solvable we must have the physical condition: 

(2.14) Vat/i<-c o <0, on dQ, where V A f=M i d xt . 

We assume that e(h) is given smooth function of h that satisfies (2.12). The condition on the metric 
is true initially since f2 is diffeomorphic to some set in R n but we have to assume that the condition 
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(2.14) is true initially. By continuity, the conditions on the inverse of the metric and on the enthalpy 
are then true also for small times, with c\ replaced by 2ci and Co replaced by cq/2. In the iteration we 
construct we have to make sure that the iterates are small enough and the time is small enough that 
these conditions remain small for all the iterates. This will be discussed in section 11. In order to solve 
the wave equation on a bounded domain one also needs compatibility conditions on initial data. These 
conditions are so that the initial conditions are compatible with the boundary condition h\ QQ = 0. 
These conditions are: 

(2.15) D*h\ m = 0, when t = 0, for fc = 0, ...,m - 1 

For k = 0, 1 this is simply conditions on initial data and for k > 2 one can use (2.8) and (2.7) to express 
it in terms of lower time derivatives of h, so it can be calculated in terms of the initial conditions. 

(2.15) is called the m th order compatibility condition. We will assume that our initial data are smooth 
and satisfies the m th order condition for all m. This will be used to construct an approximate solution 
that satisfy the equation to all orders as t — > 0, and in fact the initial conditions will be turned into 
an inhomogeneous term that vanishes to all orders as t — > 0. This will also be discussed in following 
sections. We prove the following theorem: 

Theorem 2.1. Suppose that initial data (vo,fo,po) in (2.10)-(2.11) are smooth and the compatibility 
conditions (2.15) hold for all orders m. Suppose also that (2.13) and (2.14) hold when t = 0. Then 
there is T > such that (2.7)-(2.8) have a smooth solution x £ C°°([0,T] x Q), if e(h) is a smooth 
strictly increasing function with e(0) = 0. Furthermore (2.13) and (2.14) hold for < t < T with Co 
replaced by Co/2 and c\ replaced by 2c\. 

We remark that if x € C°°([0,T] x $7) and the conditions in the theorem hold then (2.8) has a 
solution h G C°°([0,T] x f2). Theorem 1.1 follows from Theorem 2.1 since initially we assumed that 
T>o is diffeomorphic to the unit ball, so (2.13) holds initially. We will first prove Theorem 2.1 when 
e(h) = ch is a linear functions of h. This corresponds to p = p(p) = co(p — po), where c, po > are 
positive constants. The result is true in general when the pressure is any strictly increasing smooth 
function of the density. Only the estimates for the enthalpy in terms of the coordinates have to be 
modified and we show how to do this in section 17. The reason we first pick a linear function is that in 
this case we get a linear wave equation for the enthalpy and that simplifies the estimates and makes it 
more similar to the incompressible case where we have a linear elliptic equation for the pressure. 

Let us now, define the Euler map, that will be used to find the solution of Euler's equations. A 
solution of Euler's equations is given by 0(x, h) = 0, where O = (©o, • • • , n ) is given by 

(2.16) &i(x,h) = D t Xi + dih, i = l,...,n, 
and 

(2.17) @ (x, h) = D 2 t e{h) -Ah- {diV k )d k V\ h\ m = 0. 

We have assumed that our initial conditions satisfy compatibility to all orders, i.e. there are smooth 
functions (xo,/io) satisfying the initial conditions (2.10)-(2.11) and 

(2.18) £> t fc 0(x o ,M| t=o = O, for all k > 0, h o\ m = 0- 

In the process of solving Q(x, h) = we will only consider functions (x, h) which has the same time 
derivatives as (xo, Hq) when t = and which satisfy h\ dQ = 0. Let us therefore introduce the notation 

(2.19) C^([0,T] xU) = {ue C°°([0,T] x H); D k t u\ t=Q = 0, for all k > 0} 
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and for short, let Cffi = Cffi([0,T\ x Q) and C°° = C°°([0,T] x Q). (2.18) then says that &{x ,h ) G 
Cqq and if x — xo G then it follows that also @ (x,h ) G Cgg. Then by [HI] we can solve 
the equation O (x,/i) = 0, with initial conditions (2.10)-(2.11) and boundary conditions h\ Qn = 0. 
The result in [HI] is formulated for vanishing initial conditions and instead an inhomogeneous term 
that vanishes to all orders as t — > 0. However, we can turn the problem into this by considering 
&o(h) = @o(x, h + h ) — ©o(x, ho) = — @o(x, h ), with vanishing initial data for h. This gives a solution 
in h G Cqq satisfying the boundary condition h\ &n = 0. 

Let (xo,/io) be the formal solution given above. For x satisfying x — x G Cqq, we now the define 
the Euler map <E> = (<J?i, . . . , <!>„) to be 

(2.20) $i(x) = D 2 tXi + dih, i = l,...,n 
where h = ^(x) is given implicitly by solving 

(2.21) D 2 t e(h) -Ah- {diV^V*) = 0, h\ dn = 
with initial conditions (2.11). A solution of Euler's equations is given by 

(2.22) $(x) = 

By the preceding argument we can, in fact, find a solution h to (2.21) such that h — h £ if 
x — xq 6 C^. The reason we choose to consider the map $(x) instead of Q(x,h) is that we must 
make sure that h\ Qn = and that the physical condition is satisfied, since the linearized operator is 
not invertible otherwise. Alternatively, one could also have tried to only consider h satisfying these 
conditions, but it seems much more difficult to preserve these conditions in the smoothing process used 
in the Nash-Moser iteration. The main work will now be to prove tame estimates for the inverse of the 
linearized operator. 

In the Nash-Moser iteration we will in fact solve for 

(2.23) $(n) = $(u + x )-$(x ) 

Let F = $(x ), when t > and F = when t < 0, and for 5 > let F 5 (t,y) = F(t - S,y). Then 
F 5 G C^. We will solve for 

(2.24) $(u) = F s - F 

The Nash-Moser theorem says that if the linearized operator is invertible and we have tame estimates 
for its inverse and for the second variation of the operator then in fact we have a solution of (2.24) if 
the right hand side is small in C^. But the right hand side of (2.24) tends to zero in C°° when 5^0 
so (2.24) has a solution for some S > and hence 

(2.25) $(« + x ) = 0, 0<t<5 

As pointed out above at each step of the iteration we will only have functions u that vanish to all orders 
as t — > 0. This condition can in fact be preserved by smoothing operators. 

In order to solve 3>(x) = we must show that the linearized operator is invertible. Let us therefore 
calculate the linearized equations. If 5 is a variation in the Lagrangian coordinates, i.e. derivative w.r.t. 
a parameter when t and y are fixed. Since [S, d/dy a ] = it follows that 

(2.26) [6A] = (s^)^ = -(dM)d h 

where we used the formula for the derivative of the inverse of a matrix 5A~ l = — A~ 1 (5A)A~ 1 . It 
follows that [5- 5x k d k ,di]=0 and by (2.20) 

(2.27) D 2 t 5 Xi - (d k D 2 Xi)5x k - di(5x k d k h - Sh) = &(x)6xi - 5x k d k $ l 
It follows from this and (2.21) that we have: 
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Lemma 2.2. Let x = x(r, t, y) be a smooth function of (r, t, y) G K = [—£, e] x [0, T] x Q, e > 0, such 
that x\ r _ Q = x. Then is a smooth function of (r,t,y) G K, such that d$(x)/dr\ Q = <&'(x)8x, 

where Sx = dx/dr\ Q and the linear map L = 3>'(x) is given by 

(2.28) &(x)6xi = D 2 5 Xi + {d l d k h)5x k - d l {5x k d k h - 5h) 

where 



(2.29) D 2 (e'{h)5h) - 5x k d k D 2 t e{h) - A(5h - 5x k d k h) - 2{d i V k )d k {5V i - 5x l diV l ) = 0, Sh 



= 

da 



In order to use the Nash-Moser iteration scheme to obtain a solution of (2.13) we must show that 
linearized operator is invertible and that the inverse satisfies tame estimates: 

Theorem 2.3. Suppose that x ,h G C°°([0,T] x TT) is a formal solution at t = 0, i.e. (2.18) hold. 
Suppose that (2.13) and (2.14) hold when t = 0. Let 

(2.30) IIMIkoo = sup sup V \d"u(t,y)\ 

o<t<T y en rrz 

Then there is a T = T(xq, h ) > 0, depending only on upper bounds for |||xo||| ro +4,oo + |||^o|||ro+4,oo; 
where r = [n/2] + 1, Cq 1 and c\, such that the following hold. If x — x G Cqq, h is the solution to 
(2.21) with h-h £ and 

(2-31) T<T , \\\x - a? ||L+4,cx> < 1, 

then (2.13) and (2.14) hold for < t < T with Co replaced by co/2 and c\ replaced by 2c\. Furthermore, 
linearized equations 

(2.32) &(x)6x = 6&, in [0,T]xfl, 

where G Cq§ has a solution Sx G Cgg. T/ie solution satisfies the estimates 

(2.33) |||<fo|||a,oo < C a (|||(5$||| a+T . 0iOO + |||(5$|||o,oo \\\x - a;o|||a+2r +4,oo), a > 

where C a = C a (xq , j Cq 1 , ci ) is bounded when a is bounded. 

Furthermore <3? is twice differentiable and the second derivative satisfies the estimates 

(2.34) \\\<f>"{x)(5x,ex)\\\ a!00 

^ ^a^|||^|||a+2r n +4,oo|||^|||o,oo + |||^|||o,oo|||^|||a+2r +4,oo + |||^ — ^0 || |a+3r +6,oo || |^ x || |o,oo || | eaJ || |o,oo) 

Theorem 2.1 follows from Theorem 2.3 (or rather a version with Holder norms.) and the Nash-Moser 
theorem; Theorem 15.1. Theorem 2.3 follows from Lemma 11.3 and Theorem 12.2. More precisely, we 
will first prove an L 2 estimate 

(2.35) \\Sx(t)\\ r + \\6x(t)\\ r < C r V|||x||| r+ro+ 4- s ,oo I \\5*{T)\\ a dT, 

s=l J 



where 

(2-36) K*)||r= E ll$>(V)IU'(n) 

|a|<r 

and we will first prove this estimate for a lower order modification, L\ , of the linearized operator to 
be described below. Furthermore in Theorem 9.1 we will first prove the estimate for L 1 expressed in 
the Lagrangian coordinates, i.e. when also the vector field is expressed in the Lagrangian frame to be 
described below. 

It will follow from estimating the solution of the wave equation (2.29) that Sh will have the same 
regularity as Sx. Since 

(2.37) [d^D^Sx* = 2{d i V k )d k {SV i - Sx 1 ^) + {d k D t V i )d i Sx k + 5x l 8^(8^)8^) 
we get by taking the divergence of (2.28): 

(2.38) D 2 t dbr5x - 5x l 8i(divD t V- {8^)8^) - A(5x k 8 k h - Sh) + 2{8 % V k )d k (5V i - Sx 1 ^) 

= div(&(x)Sx) - {8i5x l )8i& - 6x l didiv<f> 

Hence if we add (2.29) and (2.38) we get 

(2.39) div(&(x)5x) = D 2 t {d\v5x + e'(h)Sh) + (diSx 1 )^ 
There is also a similar identity for the curl, see [L3] ; 

(2.40) curl($'(x)fe) = C Dt curl (D t 5x - 5x k 8v k ) + (8 i 5x k )8 J ^ k - (a i fe fc )9 l $ fc 
where Co t is the space time Lie derivative with respect to the vector field D t = (1, V): 

(2.41) C Dt <nj = D t <Tij + {dy)a l3 + (djV l )a u 

restricted to the space components. It can be integrated along characteristics since it is invariant under 
changes of coordinates: 

(2.42) D t (aialaij) =aia J b £ Dt aij, where a i a = 8x l /8y a 

We now want to modify the linearized operator by adding a lower order term so as to remove the 
last term on the right in (2.39) or else replace it by a lower order term proportional to d\v5x and Sh 
so that the equation div(<E>'(x)&c) = gives an estimate for divfe in terms of Sh which we have better 
control of than 8Sx. In [L3] we used the modification 

(2.43) Ltfx* = L 5x i - 6x l di& + <5xMiv$ 

= D 2 Sxi - {8 k D 2 Xi)Sx k - 8i(Sx k 8 k h - Sh) + fcrMivS 

where 



(2.44) 



div$ = D t dwV + Ah + (diV^djV* = D t divV + D 2 t e{h) 
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Then we get 

(2.45) div(Li<fcr) = D 2 (div5x + e'(h)Sh) + (divSx) div$ 

L\ is a lower order modification of the linearized operator, that reduces to the linearized operator at a 
solution of <&(x) = 0. We will first prove that Li is invertible using that it has a nice equation for the 
divergence. Then once we inverted Li we will obtain estimates for L\ that are so good that they can 
be used to iterate and deal with lower order modifications like Lq. However, as pointed out in [Ha] one 
do not need to invert the operator exactly but only do so up to a quadratic error. But this operator 
L or some modification of it is likely to show up at other places so it seems important to show that a 
more general class of operators are invertible and to have good estimates for them. 

Things are somewhat easier to see if we express the vector field in the Lagrangian frame: 



(2.46) 




Then 




(2.47) 


D t W a = d g V l(5V k 5x l d l V k ), 


(2.48) 


D 2 t W a = ^(D*5xi - 5x k d k D 2 tXi - 2(5V k -5x l d l V k )d k v i ) 


Since also 




(2.49) 


Dt9ab = dy" dy b ( diVj + 9jVi ^ 



the modified linearized operator (2.43) in the Lagrangian frame become 

(2.50) L{W a = g ab D 2 W b - d a ((d c h)W c - Sh) + (D t9ac - u ac ) D t W c + div$ W a . 



Also, the divergence is invariant, 

(2.51) divfe = divW = K _1 9 (AeW°) 

Now, D t does not commute with taking the divergence, if K~ x D t n = divV / 0, so we will replace it by 
a modified time derivative that does: 

(2.52) D t = D t +divV, ie D t W a = D t W a + {d\vV)W a = n' 1 D t {nW a ) 
It then follows that 

(2.53) D t divW = div D t W 

if a = In k then & = D t a = divV, see [LI]. With this notation we have D 2 = (D t +divV)(D t +divV) = 
D 2 + 2cj A + a 2 + a = D 2 + 2&D t + a - a 2 so 



(2.54) 



D 2 = D 2 - 2a D t + & 2 -a, D t = D t -& 
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Hence, with W = D t W and W = D 2 W, we can write the equation (2.50) as 

(2.55) L{W a = W a - g ab d b {(d c h)W c - Sh) - B 1 W a - B W a 
where 

(2.56) B x W a = -g ab {D t g bc - oj bc )W c + 2&W a , 

(2.57) B W a = g ab (D t g bc - to bc )&W c - (D 2 e(h) + a 2 )W a , 

Taking the divergence of (2.55) gives 

(2.58) div^W = D 2 t d\vW - A((d c h)W c - Sh) - divB^ - dwB W 
On the other hand by (2.29) 

(2.59) D 2 (e'(h)5h) - A (Sh - (d c h)W c ) + divB^ - 2&divW + div B W + (D 2 e(h) + a 2 ) divW = 
If we add (2.59) to (2.58) we get 

(2.60) div L{W = D 2 (divW + e'(h)5h) + d\v>5>d\vW 

as it should be, by (2.45). With tp = divW + e'(h)Sh we can hence alternatively write (2.59): 

(2.61) D 2 (e\h)5h)-A(5h-(d c h)W c )+dwB 1 W+dwB W-2&ip+(D 2 e(h)+d 2 )ip-dw^e'(h)5h = 



or 

(2.62) D 2 (e'(h)5h) - ASh + 5x k )d k h) + (tfSx^didkh - 2{d i V k )d k 5V i = 0, Sh 



= 

an 



We now, also want to express Lq = &'(x) is these coordinates. In order to do this we must transform 
the term 5x k d k $ l in (2.25) to the Lagrangian frame. If <I> a = & l dy a /dx l , then (Sx k d k <& % )dy a / dx l = 
W c V c $ a , where V c is covariant differentiation, see e.g. [CL]. Hence by (2.43) 

(2.63) L W a = L\W a - B 3 W a , where B 3 W a = -W c V c <S> a + W a div$ 

3. The projection and the normal operator. 
The energy and curl estimates. 

Let us now also define the projection P onto divergence free vector fields by 
(3.1) PU a = U a -g ab d bPu , A Pu = divU, Pu\ dn = 

(Here Aq = K~ 1 d a (Kg ab d b q) . ) P is the orthogonal projection in the inner product 



(3.2) (U,W)=f 



g ab U a W b Kdy 



a 
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and its operator norm is one: 
(3.3) 



For a function / that vanishes on the boundary define AfW a = g ab AfW b , where 

(3.4) A f W a = -d a ((d c f)W c -q), A((d c f)W c -q)=0, q\ m = 

This is defined for general vector fields but it is symmetric in the divergence free class. If U and W are 
divergence free then 

(3.5) (U,A f W)= f n a U a (-d c f)W c dS 

Jdn 

where n is the unit conormal. If f\ gn = then —d c f\ dn = (— Vjv/)n c . It follows that Af is a symmetric 
operator on divergence free vector fields, and in particular 

(3.6) A = A h 

is positive since we assumed that — V/y/i > c> on the boundary. We have 

(3.7) \{U,A f W)\ < \\V N f/V N p\\ L o am) (U,AU) 1 / 2 (W,AW) 1 / 2 

Note also that Af only depends on V/v/ on the boundary so we can replace / by something with 
the same first order derivative that is supported in a neighborhood of the boundary. We now want 
to estimate the norm of Af. Now the projection has norm one so we can drop q in (3.6). If S is a 
tangential vector field then S a d a ((d c f)W c ) = (d c Sf)W c + (d c f)CsW c , where the Lie derivative, Cs 
is defined in the next section. Furthermore if R is the normal vector fields then we can replace / by 
the distance d to the boundary times the value of Vjv/ at the boundary extended to be constant along 
the normal. Replacing / by this function we see that d c f is then independent of the radial variable so 
R a d a ((d c f)W c ) = (d c f)RW c . In conclusion we get that 

(3.8) \\A f W\\ < CJ2 II VjvSyiUcopn) || W\\ + C\\V N f\\ La o (8n) \\dW\\ 

where S is a spanning set of tangential vector fields. 

Let us now also define the projected multiplication operators Mp with a two form (3 by 

(3.9) M p W a = P((3 ab W b ) 
Since the projection has norm one it follows that 

(3-10) \\MpW\\ < \\P\U\W\\ 

Furthermore we define the operator taking vector fields to one forms 

(3.11) GW a = M g W a = P(g ab W b ) 

Then G acting on divergence free vector fields is just the identity /. 
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Let Li be the modified linearized operator in (2.55). We now want to project the equation 

(3.12) L 1 W = F 

to the divergence free vector fields: We will decompose L\ into 3 parts. We write 

(3.13) W = W + W U W = PW, W? = g ab d bqi , 

Then if g ab = D t g ab , where D t = D t - a, we have d a D t qi = D t (g ab W b ) = g a bW\ + g ab W\ and 
d a D\q x = g ab W\ + 2g ab W b + g ab W\. Hence 

(3.14) W? = g ab d b D 2 t qi - 2g ab g bc W 1 c - g ab g bc W? 

Since the projection of a gradient of a function that vanishes on the boundary vanishes it follows that 

(3.15) PW? = B 2 (W 1 ,W 1 ) a , where B 2 {W u W 1 ) a = - P (2g ab g bc Wf + g ab g bc W?) 

Let Lio = PL\. Since divWo = and Sh vanishes the boundary it follows from projecting (2.55) that, 

(3.16) L W W = W + AW - B W W - B 00 W 

(3.17) L 10 W 1 = AW 1 - BuW! - B 01 W! 

where 

(3.18) B l0 W= PB t W, B u W a = PB 1 W a + 2P(g ab g bc W c ) B 01 W a = PB W a + P{g ab g bc W c ). 
Hence the projection of (3.12) becomes 

(3.19) L 10 W = -L 10 W 1 + PF, 
Here, by (2.60) 

(3.20) W? = g ab d b q u A qi = v - e'(h)Sh, Qi\ dn = 
where 

(3.21) ti*<p + div $ (p = div F + div $ e' (h) Sh, 
and by (2.59) 

(3.22) D 2 (e'(h)Sh) - ASh = A((d c h)W c ) - div 5^ + 2&dwW - dwB W - (D 2 e(h) + & 2 ) divW 

Hence we have obtained a system of equations for (Wo, W\, (p, Sh). 

We now want to show how to show the main idea of how to obtain estimates for the divergence free 
equation. Let W = D t W = D t W + (divV)W = k^D^kW) = D t W + &W Then W is divergence free 
if W is divergence free. Let us now derive the basic energy estimate for equations of the form 

(3.23) W + AW = H 
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where A is either the normal operator or the smoothed out normal operator. Now, for any symmetric 
operator B we have 



(3.24) — (W,BW) = 2(W,BW) + (W,BW) 

where B is the time derivative of the operator B considered as an operator from the divergence free 
vector fields to the one forms, see section 4. B is defined by (4.24) with T = D t . For the two operators 
that we will consider here this is given by (4.39) and (4.40). Note that, the projection in (4.24) comes 
up here since we take the inner product with a divergence free vector field in (3.24). Let the lowest 
order energy Eq = E(W) be defined by 

(3.25) E(W) = (W,W) + (W,(A + I)W). 
Since (W, W) = (W, GW), it follows that 

(3.26) E = 2(W, W + (A + I)W) + (W, GW) + (W, (A + G)W) 



It follows from (4.24)-(4.36) and (3.7) and (3.10) that 

(3.27) |<WUW0| < Wh/hWocWAW), \(W,GW}\ < ||^||oo<w;w) 
The last two terms are hence bounded by a constant times the energy so it follows that 

(3.28) \E Q \ < V^(2||#|| +cV^), c = || V^Hoo + H^Hoo + 2 



Now, let w = W_, w = W and w = W_, i.e. w a = g a bW b , w a = g a bW b and w a = g a bW b . Observe 
that w ^ D t w etc. since W = D t W . In fact 

(3.29) D t w a = g ab W b + w a , D t w a = g ab W b + w a 
where g ab = Dtg a b and W a = D t W a . Now our equation says that 

(3.30) w + AW = H, H = —B W — B\W + F 
Since curlAVF = it follows that 

(3.32) | cuxlti>| < C(\dW\ + \W\ + \dW\ + \W\ + | curlF|) 

and hence 

(3.32) lAcurlu-l < C(|curlu;| + \dW\ + \W\), 

(3.33) |Acurlii;| < C(\dW\ + \W\ + \dW\ + \W\ + |curlF|) 
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4. The tangential vector fields and Lie derivatives 



Following [LI], we now construct the tangential vector fields, that are time independent expressed in 
the Lagrangian coordinates, i.e. that commute with D t . This means that in the Lagrangian coordinates 
they are of the form S a (y)d/dy a . Furthermore, they will satisfy, 

(4.1) d a s a = 0, 

Since Q is the unit ball in R ra the vector fields can be explicitly given. The vector fields 

(4.2) y a d/dy b - y b d/dy a 

corresponding to rotations, span the tangent space of the boundary and are divergence free in the 
interior. Furthermore they span the tangent space of the level sets of the distance function from the 
boundary in the Lagrangian coordinates 

(4.3) d(y) = 6]st(y,dO) = l-\y\ 

away from the origin y ^ 0. We will denote this set of vector fields by So We also construct a set 
of divergence free vector fields that span the full tangent space at distance d(y) > d and that are 
compactly supported in the interior at a fixed distance do/2 from the boundary. The basic one is 

(4.4) h{y\ y n ) (/ (y 1 )*/ (y 2 ) W - /'(y W) d/dy 2 ) , 

which satisfies (4.1). Furthermore we can choose f,g, h such that it is equal to d/dy 1 when \y % \ < 1/4, 
for i = 1, ...,n and so that it is when \y l \ > 1/2 for some i. In fact let / and g be smooth functions 
such that f(s) = 1 when |s| < 1/4 and f(s) = when |s| > 1/2 and g'(s) = 1 when \s\ < 1/4 and 
g(s) = when |s| > 1/2. Finally let h(y 3 , ...,y n ) = f(y 3 ) ■ ■ ■ f{y n )- By scaling, translation and rotation 
of these vector fields we can obviously construct a finite set of vector fields that span the tangent space 
when d> do and are compactly supported in the set where d > do/2. We will denote this set of vector 
fields by S\. Let S = So U S\ denote the family of tangential space vector fields and let T = S U {D t } 
denote the family of space time tangential vector fields. 

Let the radial vector field be 

(4.5) R = y a d/dy a . 
Now, 

(4.6) d a R a = n 

is not but for our purposes it suffices that it is constant. Let 7l = 5U {R}. Note that 1Z span the full 
tangent space of the space everywhere. Let U = S LI {R} U {D t } denote the family of all vector fields. 
Note also that the radial vector field commutes with the rotations; 

(4.7) [R, S]=0, S £ So 

Furthermore, the commutators of two vector fields in So is just ± another vector field in So- Therefore, 
for i = 0, 1, let TZi=SiU {R}, % = SiU {D t } and Ui = S { U {R} U {D t }. 
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Let U = {Ui}f£ 1 be some labeling of our family of vector fields. We will also use multindices 
/ = (ii, ...,i r ) of length |/| = r. so U 1 = • • • Ui r and C\j = Cjj ii • • • Cjj ir , where Cjj is the Lie 
derivative, defined below. Sometimes we will write L\j, where U G So or / G So, meaning that Ui k G So 
for all of the indices in /. 

Let us now introduce the Lie derivative of the vector field W with respect to the vector field T; 
(4.8) C T W a = TW a - {d c T a )W c 

We will only deal with Lie derivatives with respect to the vector fields T constructed in the previous 
section. For those vector fields T we have 



(4.9) [D t ,T], and [D t ,C T ]=0 
The Lie derivative of a one form is defined by 

(4.10) C T a a = Ta a + (d a T c )a c , 

The Lie derivative also commute with exterior differentiation, [Cx,d\ = so 

(4.11) C T d a q = d a Tq 

if q is a function. The Lie derivative of a two form is given by 

(4.12) C T p ab = Tp ab + (8 a T c )p cb + (d b T c )p ac 

Furthermore if w is a one form and curlu;^ = dw ab = d a w b — d b w a then since the Lie derivative 
commutes with exterior differentiation: 



(4.13) Ct cuv\w ab = cui\CTW ab 
We will also use that the Lie derivative satisfies Leibniz rule, e.g. 

(4.14) C T {a c W c ) = (C T a c )W c + a c C T W c , C T {(3 ac W c ) = (C T (3 ac )W c + (3 ac C T W c . 
Furthermore, we will also treat D t as if it was a Lie derivative and set 

(4.15) C Dt = D t 

Now of course this is not a space Lie derivative. It can however be interpreted as a space time Lie 
derivative. But the important thing is that it satisfies all the properties of the other Lie derivatives we 
are considering. The reason we want to call it Co t is simply a matter of that we will apply products 
of Lie derivatives and D t applied to the equation and since they behave in exactly the same way it is 
more efficient to have one notation for them. 

In [LI] we used extensively that the Lie derivatives with respect to the vector fields above preserved 
the divergence free condition. This is no longer true if k is not a constant, since divll = K 1 ~ 1 d a (KU a ). 
This is no longer the case if U is not divergence free. One could modify the vector fields by multiplying 
them by re -1 . However, instead we will essentially multiply the vector field we apply them to with n. 
The modified Lie derivative is now for any of our tangential vector fields defined by 



(4.16) 



C V W = C V W + (divU)W, 
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They preserves the divergence free condition, in fact 



(4.17) divCuW = UdWW, where Uf = Uf + (divf/)// 

if / is a function. This definition is invariant and (4.17) holds for any vector field U. However, in 
general, since we are considering Lie derivatives only with respect to the vector fields constructed above 
and only expressed in the Lagrangian coordinates it is simpler to us the definition 

(4.18) C V W = k~ 1 C v {kW) = C V W + (Ua)W, where a = In k 

Due to (4.1), divS = Sa if S is any of the tangential vector fields and divi? = Ra + n, if R is the radial 
vector field. For any of the tangential vector fields it then follows that 

(4.19) dwCuW = UdivW, where Uf = Uf+ (Ua)f. 

This has several advantages. The commutators satisfy \L\j,Lt\ = £[u,t], since this is true for the usual 
Lie derivative. Furthermore, this definition is constant with our previous definition of D t . 

However, when applied to one forms we want to use the regular definition of the Lie derivative. 
Also, when applied to two forms most of the time we use the regular definition: However, when applied 
to two forms it turns out to be sometimes convenient to use the opposite modification: 

(4.20) £ T (3 ab = C T p ab - (Ua)(3 ab , U = U- (Ua) 

We will most of the time apply the Lie derivative to products of the form a a = f3 ab W b : 

(4.21) C T {(3 ab W b ) = (C T (3 ab )W b + p ab C T W 

since the usual Lie derivative satisfies Leibniz rule. Using the modified Lie derivative we indicated 
in [L2] how to extend the existence theorem in [LI] to the case when k is no longer constant, i.e. 
D t a = divV / 0. This will be carried out in more detail here. 

We will now calculate commutator between Lie derivatives and the operator defined in the previous 
section, i.e. the normal operator and the multiplication operators. It is easier to calculate the commu- 
tator with Lie derivatives of these operators considered as operators with values in the one forms. The 
one form w corresponding to the vector fields W is given by lowering the indices 

(4.22) w a = W a = g ab W b 

For an operator B on vector fields we denote the corresponding operator with values in the one forms 
by .B. These are related by 

(4.23) BWa = 9abBW b , BW a = g ab B a 

Most operators that we consider will map onto the divergence free vector fields so we will project the 
result afterwards to stay in this class. Furthermore, in order to preserve the divergence free condition 
we will use the modified Lie derivative. If the modified Lie derivative is applied to a divergence free 
vector field then the result is divergence free so projecting after commuting does not change the result. 
As pointed out above, for our operators it is easier to commute Lie derivatives with the corresponding 
operators from the divergence free vector fields to the one forms. Let Bt be defined by 

(4.24) B T W a = P(g ab {C T BW b - B b C T W)) 
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In particular if B is a projected multiplication operator B_ a W = P((3 ab W b ) = (3 ab W b — d a q, where q 
vanishes on the boundary is chosen so that div BW = then 

(4.25) C T B a W = (5 ab C T W b + (C T (i ab )W b + d a Tq 

and if we project to the divergence free vector fields then the term d a Tq vanishes since if T is a tangential 
vector field then Tq = as well. It therefore follows that Bj< is another projected multiplication 
operator: 

(4.26) B T W a = P((C T (3 ab )W b ) 

In particular, we will denote the time derivative of an operator by B = Br> t and for a projected 
multiplication operator this is 

(4.27) BW = B Dt W = P{{D t p ab )W b ) 

If B maps on to the divergence free vector fields 

(4.28) C T BW a = Crig^BJV) = {C T g ab )B a W + g ab C T B a W 

Here C>T9 ab = ~9 ac g bd ^'T9cd, an d if B maps onto the divergence free vector fields then £j>B is also 
divergence free so the left hand side is unchanged if we do so and we get: 

(4.29) t T BW a = -P(g ab (C T g bc )BW c ) + P{g ab {C T B a W - B^CtW)) + B£ T W a 

By (4.26) applied the G ab = P{g a bW b ) we see that G T W = P{{g ab C T g bc )W c )) so the first term in the 
right of (4.29) is GxBW a . The second term is by definition (4.24) BtW so we get 

(4.30) C T BW = BC T W + B T W - G T BW 

The most important property of the projection is that it almost commutes with Lie derivatives with 
respect to tangential vector fields: I.e. let Pu a = u a — d a pu- Then 

(4.31) PC T Pu a = PC T u a 

since CxdaPu = d a Tpjj vanishes when we project again since Tp v vanishes on the boundary. We 
have just used this fact above. We have already calculated commutators between Lie derivatives and 
the projected multiplication operators so let us now also calculate the commutator between the Lie 
derivative with respect to tangential vector fields and the normal operator. Recall that the normal 
operator is defined by AjW a = g ab Aj-W b , where 

(4.32) A f W a = -d a ((d c f)W c -q), A((d c f)W c -q)=0, q\ m = 

and / is function that vanished on the boundary. Hence since the Lie derivative commutes with exterior 
differentiation: 



(4.33) 



C T A f W a = -d a ((d c f)C T W c + (d c f f)W c + (d c Ta)fW c - Tq) 
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However, now since q vanishes on the boundary it follows that Tq also vanish on the boundary and so 
does (d c Ta)fW c . Therefore the last two terms vanish when we project again so we get 

(4.34) P(g ab jr T A f W b ) = P{g ab A f C T W b ) + P{g ab A Tf W b ) 

Let us now change notation so A = A^, where h is the enthalpy, see section 3. Then we have just 
calculated At defined by (4.24) to be At = Af h i, i.e. 

(4.35) A T = A Th , if A = A h 

In particular, if T = D t is the time derivative we will use the notation A = Ac t which then is 

(4.36) AW = A Dt W = A Dfh W 

We can now also calculate higher order commutators: 

Definition 4-1- If T is a vector fields let Bt be defined by (4.24). If T and S are two tangential vector 
fields we define Bts = {Bs)t to be the operator obtained by first using (4.24) to define B$ and then 
define (B$)t to be the operator obtained from (4.24) with B$ in place of B. Similarly if S 1 = S l2 • • • S lr 
is a product of r = |/| vector fields then we define 

(4-37) B I = (.--{B sil ).--) sih 

If B is a multiplication operator BW a = P(g ab (3 bc W c ) then 

(4.38) BjW = P{g ab (£ T Pbc)W c ). 

In particular if G is the identity operator GW a = P(g ab g bc W c ) then 

(4.39) G I W = P(g ab (C I T g bc )W c ). 
If A is the normal operator then 

(4.40) AjW a = P(g ab d b ({d c f I p)W c ) ) 

With Bt as in (4.4) we have proven that if B maps onto the divergence free vector fields then 

(4.41) C T BW = BW T + B T W - G T BW, W T = C T W 
Repeating this, gives for a product of modified Lie derivatives: 

(4.42) £ T BW = c h - IkI G h • • • G Ik B h W h Wj = C J T W 

where the sum is over all combinations of I = I\ + ... + Ik, and c Il ' IkI are some constants, such that 
C i-L...i k i — i jf j 1 _|_ j 2 — J L e ^ us then also introduce the notation 

(4.43) G hhI = c h - IkI G l3 ■■■G Ih , 

where the sum is over all combination such that Is + ...Ik = I — h — h- With this notation we can 
write (4.42) 

(4.44) C T BW = G hhI B h W l2 
where again Gj 1 i 2 j = 1 if I\ + 1% = I ■ Also let 

(4.45) G h - IkI = 0, if I 2 = I, and G h - IkI = G h - IkI , otherwise. 
Then we also have 

(4.46) t T BW = BWi + G hhI B h W h 
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5. Estimates for derivatives of a vector field in terms of the curl, 
the divergence and tangential derivatives or the normal operator. 

The first part of the lemma below says that one can get a point wise estimate of any first order 
derivative of a vector field by the curl, the divergence and derivatives that are tangential at the boundary. 
The second part say that one can also get estimates in I? with a normal derivative instead of tangential 
derivatives. The last part says that we can get the estimate for the normal derivative from the normal 
operator. The lemma is formulated in terms in the Eulerian frame, i.e. in terms the original Euclidean 
coordinates. Later we will reformulate it in the Lagrangian frame and then also get estimates for higher 
derivatives in similar fashion. 

Lemma 5.1. Let N be a vector field that is equal to the normal M at the boundary dV t and satisfies 
\N\ < 1 and \dN\ < K. Let gW = 5 ij - N i N j . Then 

(5.1) \d(3\ 2 dx < C(q kl 5 ij d k (3i dipj + | curl/3| 2 + | div/3| 2 ) 

(5.2) f \d(3\ 2 dx <C f (5 ij N k N l d t [3 k djft + | curl/?| 2 + | div/3| 2 + K 2 \f3\ 2 ) dx 
Jv t Jv t 

Suppose that S^aj is another vector field that is normal at the boundary and let Afii = di{a k (3 k — q) 
and q is chosen so that d\vA(5 = and q\gn = 0. Then 

(5.3) f 5 i3 a kai diP k dj(3 l <C f (^AfaApj + |q| 2 (| curl/3| 2 + | div/?| 2 ) + \da\ 2 \[3\ 2 ) dx 

This lemma was proven in [L3]. We will now state some results that were proven in [L3] and give 
some further definitions. 

Definition 5.1. For V and V' any of the family of vector fields introduced in section 4, and for f3 a two 
form, a one form, a function or a vector field we define 

(5-4) |/3| r v = £ \C{,0\, |/3| r v s = £ \C\jDlf3 1, 

\I\<r,IeV \I\<r,l€V,j<s 

where L\j is a product of Lie derivatives. Furthermore let 

(5.5) |/3| r = Y, l^ffl. \P\r,s= E \ d v D ^\ 

\a\ + k<r \a\<r,k<s 

where \@\ is the point wise norm. Furthermore let [P]^ = 1, [P] = 1 an d 

(5.6) [p] v B = Yl iC [p].= E 

Si + ... + S k <S,Si>l s 1 +...+s k <s,s i >l 

We note that if (3 is a function then Cjj/3 = U/3 and in general it is equal to this plus terms 
proportional to (5. Hence (5.4) is equivalent to ^|/|< r lev W 1 fl\ ■ I n particular if 1Z denotes the family 
of space vector fields then \j3\^ is equivalent to \f3\ r with a constant of equivalence independent of the 
metric. Note also that if (3 is the one form /3 a = d a q then C[j(3 = dU I q so \dq\^ = ^| 7 |< r IeV \dU I q\. 
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Definition 5.2. Let c\ be a constant such that 



(5.7) Yl (M + \9 ab \) < cl \dx/dy\ 2 + \dy/dx\ 2 < c\ 

a,b 

Let K i0 respectively Ki, for i > 1, denote continuous increasing functions of 

(5.8) ci + |||x|||i i0 ,oo, respectively || |x|| | i)00 + ci + T" 1 , 

which also depends on the order r of differentiation. Here, the norms are as in Definition 5.3. 

We note that the bounds for the metric g and its inverse follows from the bounds for the the Jacobian 
of the coordinate and its inverse since g a \, = Sij(dx l /ch/ a )(<9x J /dy b ) and g ab = 5 t3 (dy a /dx l )(dy b /dx^). 
It also follows that we have a bound for k = det (dx/dy) = y / detg and its inverse. Furthermore, if we 
have a bound for k _1 the bound for the inverse of the metric follows from the bound for the metric, 
since an upper bound for the eigenvalues and a lower bound for their product gives a lower bound for 
the eigenvalues, since they are all positive by assumption. Moreover we note that also a lower bound 
\dx/dy\ 2 > cj~ 2 follows. 

In what follows it will be convenient to consider the norms of C\jW = k _1 C\j(kW) if W is a vector 
field and of C^g = nC I u (K 1 ~ 1 g), if g is the metric. The reason for this is simply that d\v(C\jW) = 
U 1 divW and C^cuxlw = cuxl^ljW) and when we lower indices w a = g a bW b = (K~ 1 g a b)(KW b ) and 
apply the Lie derivative to the product we get Cuw a = {£u9a,b)W b + g a bCuW b . 

The following Lemma was proven in [L3] for the cases without time derivatives. The cases with 
time derivatives follows from these by applying them to W replaced by D^W, noting that dwD t W = 
D t divW and cuil(D^W)i, is equal to .D^ curly; plus terms that are lower order in time. This is similar 
to the calculation in [L1,L3] that curl^gW)^, where (CgW)^ = g a b£sW b , is equal to Cgcuilw plus 
terms of lower order. We have: 

Lemma 5.2. Let W be a vector field and let w a = g a bW b be the corresponding one form. Let k = 
det (dx/dy) = \/det g. Then 

(5.9) M + <K W , \U I K\ + \U I K- 1 \<K w c h ... Ik \U h g\---\U Ik g\ 

where the sum is over all I\ + ... + 1^ = I and K U) is as in Definition 5.2. 
With notation as in Definition 5.1 and section 4 we have 

r-l 

(5.10) \kW\? < tfio (| curH?_i + l/sdivWl^! + \kW\? + \g/ K \f_ s \ K W\f) 
We also have 

r 

(5.11) \kW\? < K 10 J2[9/^\Zs(\ cmlw ^-i + + \kW\S), 

s=0 

where the first two terms in the sum should be interpreted as if s = 0. 

The inequalities (5.10)-(5.11) also hold with (TZ,S) replaced by (14, T) and moreover, 

(5.12) \kW\?; s < K w (\ curH?_ ljS + \KdWW\?_ ljS + \kW\* s + £ \g/K\¥ +a + j \KW\?; j ) 

i+j<r+s — l,i<r,j <s 
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and 
(5.13) 



< K w £ [»/«]" +a _ i _,(|curl«;|??i 1)3 . + + \kW\fj), 



i<r,j<s 



where the sums are only over positive indices and if i = then the first two terms should be interpreted 
as 0. 

Definition 5.3. For V any of the family of vector fields in section 4 let 



(5.14) 

and let 
(5.15) 



\W\\ V r = E W^uWl 
\i\<r,iev 



\W\\ r = Y, W D td^W\\ 



\w 



\r,oo = E W C U W \ 
|/|<r,7GV 



|W|| r ,oo= E ll^^ll 
|a| + fc<r 



where ||W|| = ||W(t)|| = ||W(t, 011^(0), II^U = H^WHco = -)|U~(n)- Let the mixed norms be 

defined by 



(5.16) ||W||; 
and 

(5.17) \\W\\V ai00 

Furthermore, let 

(5.18) \\\a\" 



V 

r J 



\W\\r,, = Y E H^'^^l 
j=0 \a\<r 



Ymww^, \\w\ 

3=0 



r,s,oo = Y E 
i=0 |a|<r 



r,oo 



sup ||a(t, ' 



a 



0<t<T 



and 
(5.19) 



[ML 



ri- 



E 

.+r fc <r, 7"j>l 



r,s,oo 



ri ,oo 



sup \\a(t, ■ 



0<t<T 



r fc ,oo 



It follows from the discussion after Definition 5.1 and the beginning of Lemma 5.3 that ||W|| r is 
equivalent to ||W||^ with a constant of equivalence depending only on the dimension. As with the 
point wise estimates it will sometimes be convenient to instead use ||jCy-W|| = ||k _1 £^(kM^)||. This 
in particular true for the family of space tangential vector fields S. However instead of introducing 
a special notation we then write ||«W||f , since k is bounded from above and below by a constant 
K\ this is equivalent with a constant of equivalence K\. Furthermore, by interpolation ||reW||^ < 
^i(lllflllr||W|| + ll^llr) an d \\W\\f < Ki(\\\g\\\ r \\W\\ + ||/cW||^), and in our inequalities we will have 
lower order terms of this form anyway. Note also that it follows from the interpolation inequalities in 
Lemma 5.7 that 



(5.20) 



|[<7]jkoo< [[[0]]] PfOO <*i 



Here the constant K\ depends on a lower bound for T > and in fact tends to infinity as T — > 0. Most 
of the time this does not matter, but when estimating lower norms of the enthalpy we can not use this 
estimate. 

The following lemmas were proven in [L3] for the cases without time derivatives. Their generaliza- 
tions to including time derivatives is immediate. 
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Lemma 5.3. We have with a constant K w as in Definition 5.2, we have 

r-l 

(5.21) \\KW\\ rt0 < K 10 (\\cax\w\\ r - li0 + ||/edivW , || r _i i o + \\KW\\f fi + K w || 

.00 

s=0 

and 

r 

(5.22) \\kW\\ T)Q < K W Y^\\ S'l|r-s,o,oo(||curlu;|| s _i j o + ||/cdiv W|| s _i i0 + ||«W"||f )0 ) 

s=0 

where the two terms with s — 1 should be interpreted as when s = 0. We also have with a constant K\ 
as in Definition 5.2, we have 

r-l 

(5.23) \\kW\\ t < /fi (||curl«;|| r _i + ||KdivW|| r _i + \\kW\\J + K 10 J2 [[[9]]] r _ Si J\^W\\ s ) 

s=0 

and 

r 

(5.24) \\ K W\\ r < ^^[[[gjjj^^dlaalwW.-! + ||KdivW|| fl _i + \\kW\\J) 

s=0 

where the two terms with s — 1 should be interpreted as when s = 0. Furthermore 

(5.25) \\KW\\ r , a <ifio(||cuxl«;|| r _i, s +||«divW|| r _i >a +||KW||J a + £ [[[»]]],.+„-*-,■, JWIm) 

and 

(5.26) \\ K W\\ r , a < ^io^MU^-^ (Hcurl^H,.^ + ||«div + \\kW\\^) 

2<r, j<s 

where the two terms with i — 1 should be interpreted as ifi = 0. 

The following lemma, proven in [L3], gives a bounds of derivatives of a vector field by the curl, the 
divergence and the normal operator. We have: 

Lemma 5.4. Let cq > be a constant such that |V/v/i| > Co > 0, let K[ and K' 2 be constants such that 
||Vjv/i||too (a n) + ci < K[ and J2 S £S HViv<S7i||L°°(afi) + \dg\ +c 1 +K' 1 <K' 2 Then 

(5.27) c \\dW\\ < (\\AW\\ + K[(\\ curlu/|| + || divW||) + ^ll^ll) 

By Lemma 5.4 we also have 

(5.28) c \\dC^W\\ < ^(||cuxl(£^W)|,|| + ||div££W|| + ||AC£W|| + ||££W||) 

Here (C'j.W)^ = g a bCj,W b , and as in the proof of Lemma 5.3 in [L3] we see that curl(£j,W)|, is equal 
to curl£j,?x> = £j,curlttj, plus terms of lower order, where w a = g a bW b . In particular we see that we 
can get any space tangential derivative in this way so we also get: 
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Lemma 5.5. With K' 2 as in Lemma 5.4 we have 



r-l 



(5.29) coH^IUo < K' 2 {\\ curlu»|| r _ 1>0 + || divW|| r _i, + \\W\\f_ 1A + ||^||oo,r-»ll wr ll»,o) 

s=0 

where 

(5-30) \\W\\i A = W A ^sW\\ 

\i\= 8 ,ies 

We also have 

r-l 

(5.31) c \\W\\ r < ^(||curl«;|| r _i + ||divW|| r _i + \\W\\J_ hA + \\W\lJ_j_ + |||5|||oo,r-»ll^ll») 

s=0 

where 

(5.32) W\\l A = W At T W W 

\i\= s ,ieT 

Let us now state the interpolation inequalities that we will use. For a proof of Lemma 5.6 see e.g. 
[H1,H2] for the L°° estimate and [CL] for the L 2 estimate. 

Lemma 5.6. There are constants C r respectively C r> T depending on r respectively on r and T, such 
that if (3,W G C°°([0, T] x fi) is a two form, a function or a vector field, then 

(5-33) lll/3||U,oo<a,T|||/3||loro: /r |ll/3|||^ 
(5.34) ||W|| a < C r ||WiS~ s/r ||W||* /r 

Here \\W\\ S = ||W(t)|| s . Furthermore, for (3 G C§§([0,T] x H), see (2.19), (5.33) holds with constants 
independent of T . 

The consequence we will use is that 

Lemma 5.7. There are constants C depending on T and r, such that if f, fi, a, (3,W G C°° ([0, T] x f2) 
are functions, two forms or vector fields, then 

. fc ,oo<C|||/||fc 1 |||/||| rf00 
r,oo < C(|||a|||o,oo|||/9|||»,oo + ll|0|||o,oo||Nlr,oo) 
r-s,oo\\W\\ s < C(||| / 9||| ,oo||W r || r + |||/3||| r ,oo||W||o) 
k 

l|||si,oo - - - |||/fc|||sfc,oo < || |/l || |o, 00 - - - |||/i-l|||o,oo|||/i|||r,oo|||/i+l|||o,oo - ' ' |||/fc|||o,oo 

i=l 

where r = s\ + ... + Sk- Here \\W\\ S = ||W(t)|| s - Furthermore, for a,(3,f,fi G Cqo([0, T] x Q), see 
(2.19), (5.35)-(5.38) holds with constants independent of T . 

Proof. This follows from using Lemma 5.6 on each factor and then using the inequality A s l r B x ~ s l r < 
A + B,see[L3\. □ 
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(5.35) 


lll/llkco- 


(5.36) 


ll a II Ir- 


(5.37) 


wmw 


(5.38) 


II |/l II |si, 00' 



6. Tame estimates for the Dirichlet problem. 
In this section we will give tame estimates for the Dirichlet problem 

(6.1) Aq = f, in [0,T]xfl, q\ m = 

Let us recall the definition of the mixed norms 

(6-2) \\q\\ r ,s= E W D tdylW 

|a|<r, k<s 

In the proofs that follow we will use the interpolations 

(6-3) || | s ,oo || IffH |r,oo < ^llllffllls+r.oo) || |s+1,oo || \d\\ |r+l,oo < ^2 \\ \g\\ |s+r+l,oo 

where K\ and K 2 are as in Definition 5.2. 

Theorem 6.1. Suppose that q is a solution of the Dirichlet problem, q\ gn = and W'{ = g ab dbq then 

s 

(6.4) Wi\\t+i,. + \\dq\\r + i,s < K W £ [[M]] r+ ._ i _ il00 IIAffllij + ^10 E [[[9]]] r+s+ i-j,J W ^,J 

i<r,j<s j=0 

ifr,s > 0, and 

(6-5) ||Wi|Ui, s + |k||r+2, s <^20 £ [[[V&DtgW^+jWAqWij 

i<r, j'<s 

Moreover if P is the orthogonal projection onto divergence free vector fields and W is any vector field 
then, 

(6-6) \\PW\\ r , s < K W £ [[[9\]] r+ .-i- S ,J\W\\iJ 

i<r, j<s 

Here K\q and K20 are as in Definition 5.2, 

(6-7) [M] r>00 = E HbllU,oo • • • IIMIk.oc < ^l||b|||r,oo 

n+...+r fc <r, ri>l 

and 

(6.8) [[[(% D t g)]]] r oo = E llbllln+1,00 ■ ■ • llbllk+1,00 < -^2||blllr+i,oo 

ri + ...+r k <r, r^>l 

Remark. Note that the constant in (6.4)-(6.6) are independent of T whereas those in (6.7)-(6.8) depends 
on a lower bound for T > 0. 

First two useful lemmas: 
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Lemma 6.2. Suppose that S G S and q\ 9n = 0, and 

(6.9) £ 5 VF a = g ab d b q + F a . 
Then 

(6.10) ||£ S W|| <K 10 (||divW|| + ||F||) 
Proof. Let W s = C S W. 

(6.11) / g ab W%W b S Kdy= [ W$d a q ndy + [ W$g ab F b Kdy 

Jn Jn Jn 

If we integrate by parts in the first integral on the right, using that q vanishes on the boundary we get 

(6.12) / W§d a qKdy = - f div{W s )qndy 

Jn Jn 

However divWs = S'divH^. Then we can integrate by parts in the angular direction. S = S a d a , 
S = S + divS" so we get J n (Sf) ndy = J n d a (S a fn) dy = 0, where d a S a = 0. Hence we get 

(6.13) f W%d a qndy= f div(W) (Sq)ndy 

Jn Jn 

Here \Sq\ < \dq\ < K 10 (\W S \ + \F\) so it follows that 

(6.14) \\Ws\\ 2 <K 10 \\dwW\\(\\W s \\ + \\F\\) + K 10 \\W S \\\\F\\ 

and using that the inequality A 2 < K(A + B)B implies that A < (2K + l)B we get for some other 
constant K w : 

(6.15) ||W S || < ifiodldivWH + ||F||) □ 
Lemma 6.3. Let W a = g ab w b Then 

(6.16) D k t W a = g ab D k t w b - £ i^g ab {b k t ^g bc )blW c 



i=Q 

Proof. We have D k w b = D k (k' 1 g bc kW c ) = J2i=o (?) (^"^"Vc)) D\{kW c ), which proves the 
lemma. □ 

Proof of Theorem 6.1. To simplify notation let us denote W a = = g ab d b q. If we apply CJj. to 
w a = g a bW b we get 

(6.17) d a T I q = g ab W b + c- I ^g hab W b 2 , W I = L I T W, gi ab = t I T gab 
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and the sum is over all combinations / = I\ + 12, cJ l12 are constants such that c 7li2 = if I2 = I. We 
assume that = C J s C s Dt , where J € S and |J| = r + 1. If we write T 1 = ST K , = jCsWa-, and use 
Lemma 6.2 we get since &i\W K = f K divW = f K Aq = k" 1 T x (kA^) 

(6.18) \\WjW < K w \\f K Aq\\ + K 10 5 hl2 \\g Il \\ oo \\W l2 \\ 
or if we sum over all of them 

(6.19) \\KWf r+1 , s < K 10 \\KAq\\ r , s + K 10 £ [ML +s+ i_^ JWI^ 

i<r+l, j<s, i+j<r+s 

We now want to apply Lemma 5.3 to W a = g ab dbq. Then curlw = and divVU = Aq. 

(6.20) \\KW\\ r+1 , s <K 10 \\KAq\\ rtS + *ri||W||? +lifl + tfio E [[b]]]r+ a +i-i-i,ooll' eW llM' 

i<r+l, J<s, i+j'<r+s 

(6.19) and (6.20) therefore gives for r, s > 0, 

(6.21) IWUm <^io||KAg|| r , s + K 10 ^ [ML +s+ i-^JI KW 1k;> 

i<r+l, j'<s, i+j<r+s 

Using induction it follows that 

(6.22) \\KW\\ r+lt . < K 10 Y, [[b]]] r+ .- i - i ,coll^9ll«+^oE[[M]]r + . + i-i,ooll KW llo J - 

i<r, j'<s J=0 

(6.4) follows from this and (9 a g = g a bW b . (Replacing k by 1 just causes more terms of the same form.) 

To prove (6.5) we need estimates for the last terms on the right in (6.4). By (6.4) with r = we 
have 

s s 

(6.23) H^llo,, < K 10 J2 [MU.JIAsllo,,- + K 10 J2 [MU^JI^Ik;- 

j=0 j=0 

Since D^Aq = K~ 1 d a Df[K,g ab dbq) we have 

(6.24) HZX^gH < ||A s ^^ll + ^ioX:([[[^]]] s _ J , 00 ||^ 2 ^^|| + [[[^U^^^^ll^^ll), 
and using (6.23) it follows that 

(6.25) \\AD?q\\ < K w j2 [[[d]]] s _ 3 J\^h 3 + ^10 £ [[[d]}] s+1 _ h J\d q \\o, r 

3=0 j=0 

We have J n g ab {d a q)(dbq) ndy = — f n (Aq)q ndy and there is a constant K 10 , depending just on the 
volume of Q, i.e. J Q ndy, such that \\q\\ < ifio||Ag||, see [SY]. Hence in addition we have 



(6-26) \\q\\ + \\dq\\<K w \\Aq\\, if q\ m = 

27 



(6.26) applied to Dfq in place of q together with (6.25) gives 

(6.27) Hflgllo,, < K 10 j2 [[[am^J^Woj + K^Y: [[[9]]]. +1 _ jt Jdq\\0J, 

3=0 j=0 

where the last term should be interpreted as if s = 0, and inductively it follows that 

s 

(6-28) \\dq\\ , s < K 20 J2[[ [(dg, D t g)} ] ] s _ joo \\Aq\\ 0ij 

3=0 

It remains to prove the estimates for the projection (6.6). We have W = Wq + Wi, where Wq = PW, 
and W\ = g ab dbq where Aq = di~vW and q\ dn = 0. Proving (6.6) for r > 1 reduces to proving it for 

r = by using (6.4), since R 1 Aq = div (C^W) and replacing k by 1 just produces more terms of the 
same form . (6.6), for r = and s = follows since the projection has norm 1, ||PW|| < ||W||. Since 
the projection of g ab D^wib = g ab dbD^q vanishes we obtain from Lemma 6.3: 

k-l 

(6-29) WPD^W < K 1 '£[[[g]]] k _ i \\DiW 1 \\ 

i=0 

Since also PD^W = D^Wo we have 

(6.30) (/ - P)D^W 1 = (J - P)D k t W 
and hence since the projection has norm one 

k-l 

(6.31) ||I> t *Wi|| + \\D k t W \\ < K^DfWll+K^iMjjkJlDiW^ 

i=0 

Using induction we therefore obtain 

k 

(6-33) ||Z> t fc ^o|| + n^ii < Kl J2 [[[9]]] k _ jt JDiW\\ 

3=0 

Since as before replacing k by 1 just produces more terms of the same form this proves (6.6) also for 
r = 0. (6.7)-(6.8) follows from interpolation. □ 

7. Tame estimates for the wave equation. 

Existence of solutions for a wave equation with Dirichlet boundary conditions and initial conditions 
satisfying some compatibility conditions is well known, see e.g. [H3]. The result in [H3] is stated with 
vanishing initial condition but we will reduce to that case by subtracting of a function that solves the 
equation to all orders as t — ► 0. We consider the Cauchy problem for the wave equation on a bounded 
domain with Dirichlet boundary conditions: 

(7.1) Df(e'ip) — Aip = /, in [0,T]xfl, V| af7 = 0, 

(7.2) ^Lo=^>, A4 =0 =^i 
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Here 



(7.3) At/> = ^=0. [Vd^g ab d b i^j , 

where g ab is the inverse of the metric g ab and detg = det{<7 a {,} = k 2 , in our earlier notation. We 
assumed that e' is positive, g ab is symmetric and positive definite (since the metric is), and that g ab 
and e' are smooth satisfying: 

(7.4) e' + l/e'<c[, + \g ab \) < c?, 

a,b 

for some constants < c\ < d x < oo. 
We will apply this theorem to 



(7.5) D 2 (e'(h)8h) - A5h = -d^dx^dkh) - (PSx^didkh + 2{d l V k )d k 5V\ 5h 



= 

an 



with vanishing initial conditions and the right hand side vanishing to all orders as t — ► 0. We will 
also need some estimates for the equation for the enthalpy itself, which is best dealt with by writing 
h = h + ho where h is the smooth approximate solution. In particular, in case e(h) = c~ 2 h we get the 
equation 

(7.6) c~ 2 D 2 h -Ah = -(c- 2 D 2 h - Ah - (diV'XdjV*)) 

with vanishing initial conditions but where the right hand side vanishes to all orders as t — > and is 
smooth for t > 0. In both cases we have vanishing initial conditions. 

In the first case, (7.5), we then also want to use the special form of / in (7.1). 

(7.7) f = h + K- 1 d a (KF?) 

For the lowest order energy estimate it will be useful to take advantage of the special form (7.7) because 
it will allow us to estimate F with one less space derivative and one more time derivative instead and 
in the estimate for the divergence equation we have one more time derivative than space derivatives. 
Let 



(7.8) g s = Hlsllls+i.oo, h s = |||«||| s +2,oo, 
and let K' 2 denote a continuous function of 

(7.9) ho + go + Ci+T' 1 +r 

which in what follows also depends on the order r of differentiation. We will use that by interpolation; 

(7.10) (g s + h s )(g r + h r ) < K' 2 (g r+S + h r+s ). 
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Theorem 7.1. Suppose that initial data in (7.2) vanishes and f in (7.1) or fi and F\ in (7.7) are 
smooth and vanish to all orders as t — ► 0. Suppose also that e' = e'(h) in (7.1) is a smooth function of 
h and that h and g are smooth. Then for r > 1 the solution of (7.1) satisfies the estimates 

(7.11) |M|o,r + IMIi,r-i <i^y> r - s + <7,- s ) / (||/i||o, s -i + ||^i||o, s )dr 

a =l J ° 

Furthermore 

(7.12) U\\ r < K'^ihr-s + 9r-s){ [ (||/l || , S -1 + ||^1 ||o,») + ||/|| a _ 2 ) 

8 = 1 J ° 

and 

(7.13) U\\ r < K' 2 J> r _ s + 9r-s) l\\\f\\ S -2 + ||/||o) 

where ||/|| s -2 should be interpreted as if s = 1. 

Proof. Let us first prove (7.11). Let E r be as in Lemma 7.2 and let 

r— 1 

(7.14) El = E 2 r + f V 2 Kdy + Y, [ ga b (D s t F?)(D s t F?) ndy 

Then it follows that E r satisfy the same inequality as E r in Lemma 7.2 does, even without the last 
term in (7.19) since the norm of ip is included in E\: 

dE r 

(7.15) ~df- K '^ Y.( h rs+9r-s)(E s + WF^s + ||/l|| , S -l) 

s = l 

The highest order energy in the right hand side E s with s = r can be removed by multiplying with the 
integrating factor e Kt and integrating. We get 

r r-1 „t 

(7.16) E r < K' 2 Y^{h r _ s + g r _ s ) {\\h\\o,s-i + \\Fi\\o, s )dT + K' 2 y'(h r - s +g r - s ) E s dr 

i J o s= i Jo 

We claim that if r > 1 

(7.17) E r < K 2 J2(h r - s + g r - s ) f (||/i||o, s -i + ||Fi||o, s )dr 

If r = 1, we have just proven it in (7.16) and in general it follows from (7.16) using induction an 
interpolation. 

We claim that (7.11) follows from (7.17). In fact \\Dl~ 1 (g ab d b ^)\\ < CE r and by Lemma 6.3 
D r t -\g ab d b f) = g^d^ 1 ^ + ElZt( r - 1 )g a \Dr 1 - s g bc )D s t (g cb d b iP). If follows that ||3Z>rVll < 
Y2l=i gr- s E s which together with (7.16) and interpolation prove (7.11). 

Finally, (7.12) follows from (7.11) and Lemma 7.3. and (7.13) follows from (7.12) using that 
| s - 2 </oll/IU-2dr. 
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Lemma 7.2. Suppose that g ab and e! = e'{h) are smooth and satisfy (7.4) for r >1, 

(7.18) E r {t) = (Y,\I e '(A S+ V) 2 + 9a b (D s t ^ a )(D^ b ) K dy) 1/2 , V* = g ab d^ + F? 

s=0 ^ n 

Then for r > 1 

dE r 

(7.19) -^f^ K i H(^- s + g r -s){E s + ||Fi||o, s + ||/i||o,,-i) +K{g r _ 1 + ^_i)||V||o 

s=l 

Proof. We will prove that dE^ /(it is bounded by E r times the right hand side of (7.19), and (7.19) follows 
from this since dE r /dt = (dE 2 / dt) / \2E r ) . With the notation D t = D t + divV and D t = D t -divV we 
have D t (e'g 2 ) = (D t e')g 2 + 2e'gD t g, for any functions e' and g, and since also D t n/n = divV it follows 
that 

r-l 



JPi ft \ 

(7.20) = E/ (e / (A S+ V)(A S+2 ^)+5a b (A s ^ a )(A S+1 ^))^^ 

s=0 ^ ^ 

r — 1 

+ / ((AO(A' + V) 2 + (A^)(A** a )(A** 6 ))«dy- 

Here the terms on the second row are bounded by K' 2 E 2 . Therefore it remains to look on the terms on 
the first row. Since D t (e') = K~ 1 D t ( K Ke') it follows that 

(7.21) e'D?+^ = D?D?(e'i;)+J2B s i DW, B° = £ ^^^r^XA^'V) 

i=0 j=max(0,i-2) 

Furthermore, using Lemma 6.3 we get, since _^_ a = d a ip + F_ a , 

(7.22) Dl +1 ^ a = g ab d a D s t +1 ^ + g ab Dl +1 F b -jz{ S + l \ ab {b s t +1 ~ i g bc )b^ c 

i=o V % J 



Since e' = e'(h) it follows that the L 2 norm of the sums in (7.21) and (7.22) are bounded by 
^EUtSr-j + h r - s )E s , which is included in the right hand side of (7.19) and so is ||</ o6 .D| +1 F 6 ||. 
Therefore it remains to consider 



r-l 



(7.23) £ f ((£> t s+ V)(A S A 2 (eV) + (A'*°)(9aA a+ V))) 

r — 1 „ r — 1 „ 

= E/ ((A s+ V)(A 8 A 2 (eV)-«- 1 3„(«A 8 *°))«dj/ = E / ((A s+ V)(A s AMy 

Here we have integrated by parts using that ip\ 9n = 0, that D 2 (e'tp) — k~ x d a (n^> a } = f\ and that 
bl(K- 1 d a (n^ a )) = K- 1 d a (nb^ a ). □ 

One can get additional space regularity from taking time derivatives of the equation (7.1) and 
solving the Dirichlet problem for the Laplacian. 
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Lemma 7.3. If r > 1 then 

r 

(7.24) |M| r < K' 2 J2(9r-i + ^-i)(||^||o,i + + 11/11,-2) 

i=l 

where ||/||j-2 is to be interpreted as if i = 1. 

Proof. We can use the estimates for the Dirichlet problem: 

(7.25) = D t 2 (eV) - /, ^| 9n = 
from Theorem 6.1, since also e' = e'(/i). 

(7.26) ||^||,- SlS < ^ J] + ^-i-iJdl^IlM + ll/lk,- 2 ) 

i<r-s-2,j<s+2 

where ||/||i,j-2 is to be interpreted as if j — 2 < 0. Hence using induction and interpolation we get 

r 

(7.27) U\\ r < K> 2 J> M + /»r-i)(||V'l|0,i + Hh,i-1 + ll/IU) 

i=0 

where H^ll l is to be interpreted as if i = and ||/||i-2 is to be interpreted as if % — 2 < 0. 
Theorem 7.4. Suppose that (ft is a solution of 

(7.28) D 2 <p - A(p» = /, <f>\ t=0 = <p\ t=0 = 0. 

where p' = p'(h) = l/e'(h) and f vanishes to all orders as t — > 0; D^f\ Q = 0, for k > 0. Zei 

(7.29) W 1 = Vq, Aq = <t>, q\ m = 
and 

(7.30) F 1 =X7q', Aq' = f, q'\ m = 
Then for r > 1 

(7.31) ||W|| r+1 <^V(/i,_ s + ^_ s ) / (HFiL-i + llFilMdr 
Proof. By Theorem 7.1 

(7.32) \\<f>\\ r < K 2 J2(h r _ s + g r _ s ) f (HFilla-i + llFiydr 

»=i 

which by inverting the Laplacian in (7.29) proves that ||cWi|| r is bounded by the right hand side of 
(7.31) and it only remains to prove the estimate for ||Z)[ +1 W^i||. Using (2.54) we can write (7.28) as 

(7.33) D 2 t <j)-2d-b t (t> + k"(t>- A(p» = /, k" = a 2 - a 
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and so 

(7.34) div(Wi - V(p'dwW 1 ) - 2aW x + k"W 1 - Fi) = -2(d a &)W? + {d a k")Wf 
and hence 

(7.35) div {b 3 t Wi - D}V(p divVFi) - D{ (2aW x - fc'Wi + F r )) = -D{ (2(d a b)W? - (d a k")W?) 
We claim that 

3-1 

(7.36) Pb\W x = PB^.^Dt'W,), B,(W 1 ,...,br 1 W 1 ) = - £ $ (flp^jWi 

i=0 

In fact = PD{d a q = PD}(g ab W?) = P^Lo ® C^Sab)^!- Furthermore, let 

(7.37) Ag? = -D>(2(d a &)W? - (d a k")W?), q>\ m = 

To say that divH = is equivalent to saying that (I — P)H = so it follows from (7.35)-(7.37) that 

(7.38) b{W x = PBjiWi,..., b{- l Wx) + (/ - P)t) j ~ 2 {- V(i/divWi) + 2bW x - k"W 1 + Fi) + Vq j ~ 2 
Here 

(7.39) bt 2 (v a { P ' dw wj) = £ e7 2 )(™)(^r 2 "V 6 )^((A i -> , w)prdivw 1 )) 

and by interpolation 

(7.40) ||£r 2 (V°(j/divWi))|| <K' 2 Y, (^-2- s + /»i-3-,)||divWi|| s 

Since the projection P maps L 2 to L 2 and since inverting (7.37) maps L 2 to H 1 (in fact to i/ 2 ) it 
therefore follows that 

r r — 1 

(7.41) ||A r+1 ^ill <K' 2 Y,{9r-s \\DIW4 + {g r . 1 . s + h r _ 2 _ s ) ||divWi|| a ) + ^E^-i-sllA^ill 

s=0 s=0 

(7.31) follows from this if we use (7.32) to estimate <p = div and use that \\DIF 1 \\ < f* \\DI F x \\ dr. □ 
Corollary 7.5. With assumptions as in Theorem 7.1 we have 

(7-42) \\M\r,oo < ^(Gfc+ro-l + /»r + ro-l)lll/ll|o,oo + II l/ll |r+r -l,oo) 

where r = [ra/2] + 1. 

Proof. By Sobolev's lemma ||V'||r,oo < CllV'llr+roi which in turn can be estimated by Theorem 7.1. Then 
we estimate the L 2 norms of / by L°° norms and use interpolation. □ 

Let us now prove that the solution of (7.1) depends smoothly on parameters if the metric g and the 
inhomogeneous term / do. We have: 
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Lemma 7.6. Suppose that f G C m (B k , QjgQO, T] x 0)), G C m (B k , C°°([0, T] x 0)), w/iere 5 fe = 
{r G R fc , |r| < e}. Suppose also that g satisfies the coordinate condition uniformly in B k . Let ip be the 
solution of 

(7.43) ay = c 2 - av = /, ^| 9n = o y,| t=o = AVi =0 = o 

w/iere A is 5 wen fry (7.3). Then ip G C m (B k , QjgQO, T] x H)). 

Proof. We note that it suffices to prove the statement in the theorem for m = 1, since the general case 
follows from this by induction. In fact, assuming that V £ C' (-B fc , Coo ([0, T] x f2)), Z < m then for 
M < l: 

(7.44) QD^ = ^ c^DW^Vr, 

/3+ 7=Q , |/3|<|a| 

where 7 = (jj,...^) are ordered multi-indices with ji G {l,...,fc}. Here Z?7 = D r ij ■ ■ ■ D r -,x, where 
D r i = d/dr\ and Up ] are the repeated commutators defined inductively by \3^ l ' 1 ' > = [D r i, \3^], where 
□< ) = D r . The right hand side of (7.44) is then in C 1 (B k , Cgg([0, T] x H)) so it follows from the 
statement in the theorem for m = 1 that V G C m (B k , Cffi([0, T] x H)). 

Let us write ip r , g r f r , and d r = D 3r , to indicate the dependence of r. First we will prove that 
f r G C{B k , C^) and g r G C(5 fe , C°°) implies that ip r G C(£ fc , Cg§). We will only prove this for r = 
since the proof in general is just a notational difference from the proof for r = 0. We have 

(7.45) n r (Yv - Vo) = f r -k- (D r - Wo 

The right hand side is in C(B k ,Cqq) and tends to in C(B k ,C N ), for any iV, as r — > 0. Since in 
Corollary 7.5 we have uniform bounds for D" 1 , it follows that tp r — ipo tends to in C(B k ,C N ), for 
any N, as r -> 0, i.e. ip r G C(.B fc , Cgg). 

Let us now assume that / r G C 1 (B k , C§§) and g r G C 1 (B k , C°°). Let Vv = D r ^ r be defined by 

(7.46) D r 1p r = f r -D r 1p r , ^r-| an =0, t/V 1 1=0 = A^V 1 1=0 = 0, 

where D r = [D r ,D], f r = D r f r and D r = (D r i , D r k ). Since the right hand side of (7.46) is in 
C(B k ,C^) it follows as above that ip r G C(B k , Cqq). It remains to prove that Vv is differentiable. We 
have 

(7.47) D r ty r - Vo - W> ) =fr-fo~ rfo + - □„ - rD )^o + r(D r - Wo 

Then the right hand side divided by r tends to in C N , for any N, as r — > 0. In view of the uniform 
bounds for D" 1 in Corollary 7.5 it follows that (tp r — ipo — ripo) fr, tends to in C N for any N, as r — > 0. 
Hence V G C 1 ( J B fc ,C^). □ 

8. Tame estimates for the divergence free equation. 

Let 

(8.1) W r + AW + .BoW r + BiW' = .ff, div# = 0, VF| t=0 = W\ t=Q = 
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where H is smooth and vanishes to order r as t — ► 0, i.e. Q = 0, for k < r. Here BiW a = 

P(g ah (3 l bc W c ), i = 0, 1, are projected multiplication operators. Let 

(8.2) 58 = 11131118+1,00, ^8 = 111^1118+2,00, & = |||/3°||| a ,oo+ HI/^IIUoo, k s =g s + h s + (3 s 
and let JT^ denote a continuous function of 

(8.3) k + Cq 1 + ci + T" 1 + r, 

which in what follows also depends on the order r of differentiation. We will use that by interpolation; 
k s k r < K' 2 ' £v+s, for r, s > 0. We will prove the following estimates 

Theorem 8.1. Suppose that (2.4)-(2.5) hold for < t < T and suppose also that x is smooth for 
< t < T, T < 1. Then (8.1) has a smooth solution for < t < T. It satisfies the estimates 

r ,t 

(8.4) \\W\\ r + ||W|| r < K'i y^h-s / \\H\\ s dT, 
for r > and for m = 0, 1, 2. 

m+l 



(8.5) 



Furthermore, for r > 1, 



IIL-J- ± till 

\\D J t W\\ r < K'i E / Pt'^IU dr, 

j=0 s=0 j=0 ^° 



(8.6) ||W||r-l + ||Wl|r-l + \\W\\ r + \\W\\ r 

r—1 r t 



< K'2 ^ / (fe r _i_ a ||#|| a + A+_ s ||ij|| s + fc r+ i_ s ||F|| s + (I curlF|| s ) dr. 

s=0 J ° 



As pointed out before, existence of smooth solutions for (8.1) was proven in [L3] so we only need 
to prove the estimates. The theorem with the norms replaced by norms with just space derivatives was 
proven in [L3]. The theorem here is actually simpler to prove than the one there. 

First we note that we can reduce to the case B = B\ = since the general case follows from this 
case. Let us show this for (8.4). Assume that (8.4) holds for the case when B = B\ = 0. Then using 
(8.4) applied to the equation W + AW = Hi = H + B W + B\W gives using that, by Theorem 6.1, 
\\BiW\\ s < K^J2l =0 (f3 s _ k +g s _ k )\\W\\ k and (8.4), 

(8.7) \\W\\ r + \\W\\ r <K%J2(g r - s + h r - s + P r - s ) [ (\\H\\ s + \\W\\ s + \\W\\ s ),dT 

We can now first remove the highest order terms from the right hand side by a Gronwall type of 
argument since they occur in the left. In fact let f(t) = J* \\W\\ r + ||W^|| r <iT. Then by (8.7) /' < 
K^f + K![ H-f^llr + S r _i, where X r -i is the sum of the first r — 1 terms in the right of (8.7). Hence 
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multiplying by the integrating factor e K ' 2>t and integrating gives f(t) < K' 2 ' / * \\H\\ r + £ r -i dr. It follows 
that 

(8.8) \\W\\ r + \\W\\ r <K%y2{g r _ s + h r - s + p r - s ) [ \\H\\ s dr 

+ K'i V (<?,-. + hr-s + Pr-s) I (\\W \\ s + \\W \\ s ) , d T 
s=0 J 

This proves (8.4) for r = 0. (8.4) for r > 1 now follows by induction, using (8.4) for s < r — 1 in the 
terms on the second row of (8.8) together with the interpolation (8.4). The proof of that we can reduce 
to the case Bq = B\ = also for (8.5) follows in a the same way, using that we have already proven the 
estimate for smaller m. The prove that (8.6) can be reduced to the case B = Bi = 0, we also estimate 
the lower order terms using (8.5). 

We now in what follows assume that B = Bi = 0. Of course we could have included these operators 
in the calculations that follows but the argument becomes more clear without them. Lowering the indices 
in (8.1): 

(8.9) GW + AW=GH 

Let Clp, I £ T, stand for a product of Lie derivatives of |/| vector fields in T and let W\ = L l T W. If we 
apply repeatedly apply Lie derivatives Ct and the projection in between, see section 4, we obtain 

(8.10) c h 72 {G h W l2 + A h W l2 - G h H l2 )=0 

where the sum is over all combination of 7i + 12 = I and c /l/2 = 1. Here Gj and Aj are the operators 
given by (4.39) and (4.40). If we raise the indices again we get 

(8.11) Wi + AW 1 = -c Ill2 A Il W h - c h h G h W h + c h h G h H h 

where c hh = 1, if \I 2 \ < \I\, and c hh = if \I 2 \ = \I\. 
Let us define energies 

(8.12) E I = (W I ,W I ) + (W I ,(A + I)W I ), Ej= £ 

\i\< s ,ier 

Note that in the sum we also included all time derivatives £.D t - The reason for this is that when 
calculating commutators second order time derivatives show up in the first term on the right in (8.10). 
We get by differentiating (8.11), see the end of section 3, 

(8.13) Ej = 2(W>, Wi + (A + J)W» + <W>, GW» + <W>, (A + G)W» 

We now, want to estimate the right hand side by Ej, where r = |/|. Here, by (3.27), the last two terms 
can be bounded by (HqCq 1 +go)Ei. Therefore, (8.12) can be estimated by the L 2 norm of the right hand 
side of (8.10). In the sums with c 7li2 , in (8.9), we have I/2I < |^|- Since we included time derivatives 
up to highest orders in Ej it follows that ||W/ 2 || — . Here Gi 1 is a bounded operator so the term 
with Gi 1 can be controlled. Ai 1 is an operator of order one so the term with Ai 1 can be controlled by a 
constant times ||9Wj 2 || + || Wj 2 1| . However, at this point we only have control of tangential derivatives 
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of W. One could combine the estimate here with the curl estimate given later to get control over all 
derivatives up to order r = \I\. Instead we will add a lower order term to the energy such that the 
its time derivative cancels the terms with Aj 1 and replaces them with lower order terms that may be 
controlled. Let 

(8.14) D I = 25 I ^(W I ,A h W l2 ) 
where the sum is over all h + h = I, with|I 2 | < |/|. Then 

(8.15) Dj = 2c hl2 {(W I ,A h W l2 ) + {W!,A h W l2 ) + (W:, A h W h )) . 
Hence 

(8.16) E! + D! = 2~c^ ((Wr, A h W l2 ) + (Wj, A h W l2 }) + (W T , AWj) 

+ 2<W>, -t hl2 G h W l2 + c hl2 G h H l2 + W» + <W>, GW» + <W>, GW» 

Here, the terms on the first row can be controlled using (3.7) and the terms on the terms on the second 

row can be controlled using (3.10). We get 

(8-17) 

\(U,AjW)\ < |||/i||U + i t ooCb 1 <?/,AI7) 1 / 2 <W,AW> 1 / 2 , \(U,GjW)\ < \\\g\\\ a ,oo\\U\\\\W\\, \J\ = s 
Hence, we have proven that 

r r— 1 



'It— s c o ^ 

s=0 s=0 



r-1 



(8.18) \E I + D I \ < CE^J^^rsCo 1 + 9r-s){Ej + \\H\\J), \D T \ < CEj ^ K^E] 
so it follows that 

(8.19) Ej <C(h c^+go) [\eJ + \\H\\ r )dr + CV^-.Cq 1 + g r _ s )( [ (eJ + \\H\\J) dr + Ej) 

Jo s=0 Wo ' 

Using induction and interpolation (8.4) we get the same inequality without Ej, for s < r — 1, but 
instead multiplied by a constant K'^ depending on (8.3) and r. Then, by a Gronwall type of argument, 
see the beginning of the proof, we can also remove Ej from the right hand side of (8.17), i.e., with 
/(*) = ll E r dr, we get and inequality /' < K£f + K% Y? s =o( h r-s + 9r-s)\\H\\J and multiplying by 
the integrating factor e K ^ * gives: 

Lemma 8.2. For r > we have 

(8.20) E r r < K'2 V(/i,- s + 9r- s ) I \\Hf s dr 



This proves the estimates for tangential derivatives. We now also want to have estimate for the curl 
and then by the results in section 5 the estimates for all derivatives will follow from this. 

Let w a = g ab W b and w a = g a bW b and let cur\w ab = d a w b - d b w a . Then D t w a = w a + g a bW b , 
where g ab = D t g ab and W b = D t W b . It follows that £> t curlio a6 = cuvlw ab + d a (g bc W c ) - d b (g ac W c ) 
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Since the curl of A vanishes it follows that curlti) = curl if, if w a = g a bW b . Since the curl commutes 
with the Lie derivative it therefore follows that 



(8.21) lAcurlwli^i < C^g r - s \W\ u s + \cvx\H^_ 1 

s=0 
r 

(8.22) |Acurlu/|"_i < C^g r - a \W^ 

s=0 

and by Lemma 5.3: 



(8.23) \W\? < ^^gr-i-silcwlw^ + ldWWfJ + \W\J)+K l9r _ 1 \W\ 

s=l 
r 

(8.24) \W\ u r < K^gr-i-silcmiw^ + \ divW\^ + \W\J) + K^.^Wlo 

s=l 

Let 

(8.25) = || curlu;|| r _i + || curlu;|| r _i. 
Since divVy = divVF = it now follows that 

r r 

(8.26) \C¥\ < K^dr-siWWWs + ||W|| s ) + <7||curLH;||r-i < K x ^ g r - s {C u s + Ej) + C\\ curLff||.-i. 

s=0 s=0 

This together with Lemma 8.2 and the argument for its proof gives 
Lemma 8.3. For r > we have 

(8.27) \\W\\ r + \\W\\ r + E^ <K^J2(g r _ s + h r - s ) f \\H\\ s oIt. 

s =o J( > 

This proves the first part of Theorem 8.1. 

In order to prove the second part of Theorem 8.1 we replace the energy in (8.10) by 



(8.28) El m = V^7 

I=h+h, heT, \h\<r, I 2 = {D t ,...,D t },\I 2 \=m 

i.e. we take two additional time derivatives. Noting that, the argument leading up to Lemma 8.2 only 
requires that we have at least as many time derivatives as tangential space derivatives so it follows from 
its proof that 

Lemma 8.4. For r > we have 

r m ,.t 

(8.29) < m < K'i E kr+m-j-s / \\DiH\\J air 

s=0 j=0 ^° 
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We have w a = -A a W + H a and w a = g ab D t W b = D t w a - (D t 9ab)W b so, since curlAW = 0, 
curlu) = curlil and curlw a6 = cml(D t H_) ab - d a {b t g bc )W b ) + d b ((D t g ac )W c ) . Hence |curlu>| s _i = 
| curli7| s _i and |curlui| s _i < K\ J2 S j=o9s-j \ + \D t curliJ | s _i so using the estimates (8.23)-(8.23) 
with (W, W) replaced by (W, W) gives: 

r r 

(8.30) \W\? < K^gr-i-sW^ + K 1 Y J 9r-i-s\cm\H\ s _ 1 

S = 3 = 1 

r r 

(8.31) {WW < K^gr-i-sdW^ + \W\J) +K 1 ^2g r - 1 - 8 \D t cm\H\ 8 - 1 

s=0 s=l 

Using the estimate in Lemma 8.3, with r replaced by r — 1, the estimate in Lemma 8.4, the estimate 
\\Dt cuv\H_ || s _i < J Q * curl^ || s _i dr, and interpolation gives: 

Lemma 8.5. For r > and m = 0, 1, 2, we have 

m+l r m „( 

(8.32) ^ ll^'^llr + < m < ^2 E E i-' / II- dr - 

.7=0 s=0j=0 ^° 



We now want to use the estimate for the additional time derivatives to get estimates for and 
\\W\\ r in (8.6) since by Lemma 5.5: 

r-l 

(8.33) c \\W\\ r < ^'(Hcm-HI^! + ||divWn| r _i + ||W||^_ 1)i4 + \\W\\f-l + ^5r-l- S ||W|| s ) 

s=0 
r-l 

(8.34) c \\W\\ r < K%(\\ curlTi;)!,.! + ||divW|| r _i + \\W\\J_ hA + \\W\tf_ r + Y,9r-i-s\\W\\ s ) 

s=0 

where 

(8.35) \\W\\l A = Yl W A ^tW\\ 

\i\= s ,ieT 

The terms of order r — 1 or less in (8.33)-(8.34) can be estimated by Lemma 8.3 and Lemma 8.5 with r 
replaced by r — 1. Therefore it only remains to estimate the terms involving the curl and the operator 
A. The estimate for ||AWj|| with | J\ < r — 1 we get from (8.11); 

(8.36) AW j = -c Jl j2 A. h Wj 2 - c h J - {G. h W h - G h H. h ) 

where c Jl j2 = if | J2I = r — 1. Here the terms in the parenthesis can be estimated by Lemma 8.5, with 
r replaced by r — 1. Note that ||-ffj 2 || < J* * ||-Hj 2 || dr since we assume that H vanishes to all orders as 
t — ► 0. Since by (3.8) Aj 1 is order one and | J2I < r — 2 in the first term on the right it follows that also 
this term can be estimated by Lemma 8.3 with r replaced by r — 1. We hence get 

(8.37) \\W\\J_ 1A <K^J2 {k r -i-s\\H\\ s + k r _ s \\H\\ s + k r+1 _ s \\H\\ s ) dr 

s=0 J ° 
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From (8.9) we also get the estimate for with | J| < r — 1 by letting one of the derivatives in 

(8.9) be a time derivative so / = J + {D t }. If we write this out we get 

(8.38) AW j = -c Jl J2 A j, Wj 2 - c Jl ^ A j, Wj 2 - c Jl Ja (G Jl Wj 2 + Gj x Wj 2 - Gj, H h - Gj, Hj 2 ) 

where c JlJ2 = if | J2I = r — 1. Here the terms in the parenthesis can be estimated by Lemma 8.5, 
with r replaced by r — 1. Since Aj l and Aj 1 is order 1, see (3.8), we can also estimate all the other 
by Lemma 8.3 with r replaced by r — 1, apart from the term with AWj. This term is estimated by 
MWjII < C^o(||VF||r + ||W|| r _i), where the last term can again be estimated by Lemma 8.3 with r 
replaced by r — 1. We hence get 

(8.39) ||^||^-i a < ho\\W\\ r + K% V / (k r _ 1 _ s \\H\\ s + k r _ s \\H\\ s + k r+1 . s \\H\\ s )dT 

s=0 J 

(8.33) and (8.34) together with (8.37) and (8.39) therefore gives 

r-i ,.t 

(8.40) co||W|| r < K'2 ||curlu>|| r _i + K'^ V] / (k r - 1 - s \\H\\ s + k r - s \\H\\ s + k r+1 - s \\H\\ s ) dr 

and 
(8.41) 

co\\W\\ r < ho c o 1 II curlw|| r _i + if 2 II curlu;|| r _i +K% / (k r -i- s \\H \\ s + fc r _ s ||ii|| s + Av + i_ s ||ii|| s ) dr 

It therefore only remains to control the curl. With = || curlu;|| r _i + ||curlw|| r _i it hence follows 
from (8.26) and (8.40)-(8.41) that 

r-l ,.t 

(8.42) < K'lcy + K'l V / (Av-i-Jff|| a + *V-Jff||» + *V+i-J#y dr + || curl|i||r-i 

where we also used that, by Lemma 8.3, ||W"||o + ||W"||o < if 2 Jo H^llo^r. Integrating this equation 
gives a bound for in terms of the integral in (8.38). Using this bound in (8.40)-(8.41) then gives 

Lemma 8.6. For r > 1 we have 

r-l ,.t 

(8.43) \\W\\ r + \\W\\ r < if 2 ' V / (Av-i- s ||#|| s + AV-Jii|| s + £v+i- s ||i2"|| s + || curl^Ha) dr 

This concludes the proof of Theorem 8.1. 

9. Existence and tame estimates for the inverse of the modified linearized operator. 

We want to show existence and tame estimates for the inverse of the linearized operator. However, 
first we will show existence and estimates for the modified linearized operator L x given by (2.55): 

(9.1) L 1 W = F, W\ t=0 = W\ t=Q = 
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where F is smooth and vanishes to all orders as t — > 0. 
Let 



(9-2) n s = |||5||| s +2,oo + |||w||| s +i,oo + |||^||| s +3,oo, 

and let K% be denote a continuous function of 

(9.3) no + Cg 1 +T" 1 + ci +r 

which in what follows also depends on the order of differentiation r. In the proof that follows we will 
use the interpolation 

(9.4) n s n r < K^n s+r 

Theorem 9.1. Suppose that (2.4)-(2.5) hold for < t < T and suppose also that x is smooth for 
< t < T an that T < 1. Then (9.1), where F is smooth and vanishing to all orders as t — > 0, has a 
smooth solution for < t < T. It satisfies the estimates 

(9.5) \\W\\ r + \\W\\ r < K^J2n r _ s f \\F\\ s dT, r>l 



Proof of Theorem 9.1. The proof of existence and the estimate (9.5) for the modified linearized operator 
uses the orthogonal decomposition of the vector field into its divergence free part and a gradient of a 
function vanishing on the boundary. The solution will be further divided into four parts and for each of 
them we use either Theorem 6.1, Theorem 7.1 or Theorem 8.1. This gives us the estimates in Theorem 
9.2. These estimates imply the estimate (9.5). Furthermore, the estimates holds for iterates, i.e given an 
iterate for W we define Sh = W(x)W by (9.10) and then define W tj by (9.8)-(9.20) and (9.12)-(9.13). 
This gives us a new iterate for W. Theorem 9.2 then gives us uniform bounds for the iterates and 
applied to the equations for differences of iterates gives us convergence, see [L2,L3]. 

Now, the solution W of 

(9.6) L 1 W = F 
can be obtained as the sum of four terms, see section 3, 

(9.7) W = W + W U Wi = Wio + Wu, W = W 00 + W 01 
where Wq is divergence free and 

(9.8) W? t =g ab d b q u , i = 0, 1, 9ii| an = 0, Aq n = <p, Aq w = -e'(h)5h 
and ip satisfies an ordinary differential equation: 

(9.9) £> t V + div$<£> = dwF + div^e'{h)5h, tp\ t=Q = ip\ t=Q = 
where div<l> = D^e{h) + D t d\vV and 5h satisfies the wave equation 

(9.10) D 2 t {e'(h)5h) - ASh = A((d c h)W c ) - divB^ + 2adwW - dwB W - {D 2 t e{h) + a 2 ) divW, 

(9.11) M\ dn = 0, Sh\ t=o = Sh\ t=o = 0. 
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Here, the divergence free parts satisfy the evolution equation for the normal operator, (3.16) 

(9.12) W o + AWoo-B 10 W 00 -BooW 00 = -AW 10 + B 11 W 1 +B 01 W 1 +PF, W 00 \ t=0 = W 00 \ t=(j = 
and 

(9.13) Wot + AW 01 - B w W i - B 00 W 01 = -AW 1U Wbi| t=0 = Wbi| t=0 = 
If 

i 

(9.14) E r = WijWr + WijWr + \\5h\\ r 

i,j=0 

then we get from Theorem 9.2 that for r > 1 

(9.15) £ r <^E n " / (\\F\\s + E s )dT 

s=i J o 

Using a Gronwall type of argument and induction as in section 8 it follows that for r > 1 we have 

(9.16) E, 



f 

'r < K'l Vn r - S / \\F\\ s dT 

8=1 J ° 



which proves Theorem 9.1. □ 
Theorem 9.2. Ifr>\ we have 



(9.17) H^llllr + IIWllllr <Kg Vll r _ s / ||F|| s + ||Wi || s dT 

(9.18) ||#n|| r <K'jJ2 n r-s(\\F\\ s + \\W w \\ s + \ \\F\\ S + \\W W \\ S dr) 

r ~ r 

(9.19) \\Woi\\r + \\W 01 \\ r <K^n r _ s (\\F\\ S + \\W 10 \\ S ) dr 

8=1 J ° 

r „ r 

(9.20) ||Wbo||r + ||Woo||r <K'i Vn r _ s / (||F|| fl + ||SWi || fl )dr 

s=l J ° 

(9.21) \\W 10 \\ r + \\W 10 \\r + l <^E«" / (II^IU + II^IU)^ 

where W = W + W lf W x = W w + W n and W = W 00 + Woi . 

Proof of Theorem 9.2. (9.17) follows directly from applying Lemma 9.3 below to (9.9). Here we write 
divFi = divF + div<& e'(h)5h = div(F - div$W 10 ) + (d a div$)W? and use Theorem 6.1. This also 
gives the additional estimate (9.18). 

To prove (9.19) we use the second part of Theorem 8.1 applied to (9.13). Since curlAH^n = and 
we have an estimate for and additional time derivative in (9.18), (9.19) follows. 
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By the first part of Theorem 8.1: 
(9.22) ||Wbo||r + ||Wbo||r <^Vn r _!_ s / {\F\\ s + \\dW^\\ s + ||Wi|| s + ||Wi|| s )dr 

s=0 J ° 



Using (9.21) and (9.17) gives (9.20). Note that the sum above starts at s = but since the lower norms 
are included in the higher norms and since we replaced n r _ s _i by n r _ s we can start the sum in (9.20) 
from s = 1. 

The estimate (9.21) follows from Theorem 7.1 applied to (9.10). If / denotes the right hand side of 
(9.10) then 

S 

(9-23) ||/|| s _ 2 < K'^Y, n s-i-k{\\W\\k + \\W\\ k ) 

k=0 

Since ||/|| < K^(\\W\\ 2 + \\W\\i) we obtain from (7.13) for r > 2 

(9.24) ||<fa|| r < K'l V n r _!_ s / (|| W\\ s + \\W\\ S ) dr 

For r = 1 we use (7.12), which gives (9.24) also for r = 1 if we write / = f\ + A(h c W c ). This proves that 
||<5/i|| r is bounded by the right hand side of (9.21). To get the estimate also for Wio we use Theorem 6.1 to 
estimate the solution of the Dirichlet problem in (9.8) using that He^/i)^^ < K$ X^I=o ll^l|r-s,oc II^IU- 
This gives that ||cWio|| r and ||Wio|| r are bounded by the right hand side of (9.21) but it remains to 
show that ||L>[ +1 W r io|| is bounded by the right hand side of (9.21) which follows from Theorem 7.4. □ 

Lemma 9.3. Let ip be the solution of 

(9-25) D*ip + k<p = f, ip\ t=(j = tp\ t=Q = 0, 

where D 3 t f\\ t = for j > and k = div<E> = D^e(h) + D t divV. Then for r > 1 we have 

(9.26) N^ll, — i + IM|r-i < K'J V n r _x_ s / ||/|| s dr, 

(9-27) ll^ll^x < K'/j2 n r-i- s (\\f\\s + /||/|| a dr) 

s=0 ^° 

Furthermore let W\ and Fi be defined by 

(9-28) W 1 = Vq, Aq = ip, q\ dn = 

(9-29) F l =Vq l Aq' = f, q'\ m =0 

Then for r > 1 we have 

(9.30) ||^i||r + ||Wi||r <i^Y> r _ s twFiWsdr, 

(9.31) \\W 1 \\ r <K'^J2n r+1 . s (\\F 1 \\ s + f {{F^dr) 
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Proof of Lemma 9.3. (9.25) is an ordinary differential equation and (9.26)-(9.27) follows as in the proof 
of Proposition 10.1 in [L3]. (9.25) can be written 

(9.32) D*<p - 2&D t ip + k'ip = f, k' = a 2 + D 2 t e(h) 
and so 

(9.33) div (Wi - 2&W 1 + fc'Wi - Pi) = -2(5 a a)W 1 a + (da*/)^ 
and hence 

(9.34) div {D{W X - £>i (2&W 1 - k'W 1 + Pi)) = -P? (2(<9 a cr)W^ - (d a fc')WT) 
We claim that 

(9.35) PP^ = PB 3 (^,..,CVi), S J (^ 1 ,...,i)r 1 ^i) = E 

i=0 

In fact = PPRg = PP^T^) = P^Lo (•') G&TVab) A^i- Furthermore, let 

(9.36) Aq* = -P? (2(S 0< 7)W7 - (d a k')W") , q J \ m = 

To say that divPT = is equivalent to saying that (/ — P)H = so it follows from (9.34)-(9.36) that 

(9.37) b{W x = PBjiWi, b{~ 1 W 1 ) + (/ - P)D J t ~ 2 (2&W 1 - k'W x + Pi) + Vq j ~ 2 

Since the projection P maps L 2 to L 2 and since inverting (9.36) maps L 2 to H 1 (in fact to H 2 ) it 
therefore follows that: 

r r — 1 

(9.38) < K' ; ;Y / n r-s\\D s t W 1 \\+Y,nr-2-s\\D s t F 1 \\ 

s=0 s=0 

Sine it also follows from (9.26) that 



(9.39) 



(9.40) 



||aWi|| r _i + ||aWi|| r _i <l^Vn r _ s / ||Pi|| s dr, 

\\dW4r-! <K%J2n r - s (||Pi|| s + / ||Pi|| s dr) 



the lemma follows from also estimating \\DfFi(t, -)|| < J Q \\D^ + Pi(r, -)|| dr. □ 

10. Estimates for the enthalpy in terms of the coordinate. 

We have now proved that the linearized operator is invertible. However, since we think of h = <J>(x) 
as a functional of x we must also estimate the L°° norms of h in terms of the L°° norms of x. For 
the corresponding problem for the incompressible case in [L3] we could take advantage of the Schauder 
estimates. However for the wave equation there are no L°° estimates that do not loose regularity. 
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For wave equations it is best to get the L°° norms from the L 2 norms using Sobolev's lemma. These 
estimates where obtained in Corollary 7.5. However, in the estimates there the L°° norm also occurred 
in the right hand side, due to that we assumed that e' = e'{h). This estimate can easily be improved 
by estimating the L 2 norm of the solution of the nonlinear wave equation instead. We will however for 
simplicity assume that e(h) = ch, so that e'(h) = c, where < c < oo is a constant. In that case there 
are no h terms in the results in section 7. We have 

(10.1) D 2 ch -Ah = {diV'XdjV*), h\ dn = 

The estimates in Corollary 7.5 were however formulated for vanishing initial data. Therefore let h = 
h — ho and where ho satisfies the equation (10.1) to all orders as t — > 0, if x — xo vanishes to all orders 
as t — > 0, see section 2. It follows that 

(10.2) D 2 (ch) -Ah = -D 2 (ch ) + Ah - (diV^djV 1 ) 

vanishes to all orders as t — > 0. We can therefore apply Corollary 7.5, which gives: 
Lemma 10.1. We have 

(W.3) IHfflHa+l.oo + IIMIk oo < i^l|||x||| s+2 ,oo 

Suppose that e(h) = ch, where < c < oo. Then for r > 2 
(10-4) IH/illUoo <K 2 \\\x\\\ 

r+r +l,oo- 

Here Ki are as in Definition 5.2 and K 2 also depends on a bound for |||/io|||r+r +i,oo- 

Proof. The first inequality follows directly from the definitions and interpolation. If / denotes the right 
hand side of (10.2) then 

(10-5) IH/IIUoc < tf 2 (|||/lo||| s+2 ,co + llklll-+2,oo) 

As pointed out above if e(h) = ch then h r in section 7 vanishes and by Corollary 7.5 and interpolation 
we have 

(10-6) Hl^llkoo < -KTadN II |r+ro+l,oo II I ^0 II I r+r + l,oo) 

However since we are just looking on fixed initial data we can also include the norms of ho in the 
constants and since in fact we also have a lower bound for |||x|||i !00 > Cc\ > 0, by the coordinate 
condition, the lemma follows. □ 

Remark. Note that K 2 in (10.4) depends on ho- However, ho is a function which is fixed once we fixed 
initial data so this just leads to an r dependence of the constant. 

It now follows that with Ki as in Definition 5.2 and K[ as in sections 7, 8, and 9 we have 
(10-7) Ki < K l+ro+1 
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11. Estimates for the physical condition and coordinate condition. 

We assume that the physical condition and the coordinate condition initially at time for some 
constants cq > and c\ < oo and we need to show that this implies that they will hold with cq replaced 
by Co/2 and c\ replaced by 2ci, for < t < T, if T is sufficiently small. 

Now, for the coordinate condition this is easy since we can just estimate the physical condition by 
the time derivative of g, which can be estimated by 1 1 1 a? 1 1 1 2,00 • 

Lemma 11.1. Let M(t) — sup^gfj \J\dx / dy\^ + \dy/dx\^. Then 

(11.1) M(t) < 2M(0), for t<T, if T|| |z|| | 2jOO M(0) < 1/8 

Let N(t) = sup y€dn |V/v/i| _1 . Then assuming that T is so small that (11.1) hold we have 

(11.2) N(t)<2N(0) for t<T, if T|||/i||| 2iOO M(0)JV(0) < 1/8 

Proof. We have \D t dx/dy\ < || |a; || 1 2,00 and \D t dy/dx\ < \dy/dx\ 2 \D t dx/dy\ so M'(t) < (1 + 
M 2 )|||x||| 2! oo < 2Af 2 |||x||| 2 ,oo, since also M(t) > 1. Hence 

(11.3) M(t) < M(0)(l -2|||x|||2,ooM(0)t)" 1 , when 2|| |x|| | 2 ,ooM(0)t < 1. 

Now, V N h = N a d a h, where N is the unit normal, so D t V N h = V N D t h + (D t N a )d a h = V N D t h + 
(D t N a )g ab N b V N h, since h\ m = 0. Furthermore = D t {g ab N a N b ) = 2g ab {D t N a )N b + (D t9ab )N a N b 
and N a = (dy a /dx i )N i , where S^WNi = 1. Hence \D t V N h\ < M(\dD t h\ + \dD t x\\V N h\) Therefore 
if N(t) = sup y€dn IVat/iI -1 , we have N' < M\\ \h\\ | 2 ,oo^V 2 + M\\\x\\\ 2 ,ooN/2 and if we use (11.1) and 
multiply with the integrating factor, N(t) = N(t)e~ tM( - ^ x ^ 2 ^ we get N' < 2e 1 / 8 M(0)|||/i||| 2)OO iV 2 . 
Hence 

(11.4) N(t) < N{0)e^ s (l - iV(0)2 e 1 / 8 M(0)|||/ l ||| 2)OO t)- 1 when iV(0)2 e 1 / 8 M(0)|||/ l ||| 2)OO t < 1 
This proves the lemma. □ 

To satisfy the condition in (11.1) we just need to choose T so small that T\\ \x\\ | 2 ,ooCi < 1/8. 
Remarking that x = u + x , where xo is a fixed and that in the Nash-Moser iteration we will only apply 
our estimates to functions satisfying |||u||| r0+ 4 jOO < 1. 

However for the physical condition this is a bit more difficult. One has to control the |||/i|||2,oo 
and the estimate for h in terms of x used interpolation so they are in terms of a constant K that 
is a continuous function of T _1 and tends to infinity as T — > 0. In the compressible case we never 
used interpolation in time so the corresponding estimate there was easier. We therefore have to redo 
the estimates for the wave equation for the lowest norms without using interpolation in time. This is 
however, follows from standard estimates for the wave equation. Those we have here also work if we do 
not use interpolation. If we do not use interpolation, then the estimates in sections 5,6 and 7 still hold, 
with constants independent of T < 1, but instead of depending linearly on the highest norms of x they 
are polynomials in the highest norms of x occurring in the estimates. 

Lemma 11.2. There are continuous increasing functions C r such that for T < 1 we have 

(H-5) HHI|r,oo < C r (ci, |||a;||| T . +T . + i )OO , |||^o|||r+r +l,oo) 

Proof. The proof is the same as the proof of Lemma 10.1 using that the estimates in sections 5,6 and 
hold with constant of the form above. □ 



Summing up, we have hence proven that 
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Lemma 11.3. Let C 2 be as in Lemma 11.2 and let c\ and Co be constants such that the coordinate 
condition (2.13) and the physical condition (2.14) hold when t = 0. Suppose < T < 1 is fixed such 
that 

(11.6) Tci(l + |||so||koo)< 1/8, rCiC2(2ci,l + |||s |||ro+3,oo,|||/lo||U+3,oo) < Co/8, 

where C 2 is as in Lemma 11.2. Then for < t < T , the coordinate condition hold with c\ replaced by 
2ci and Co replaced by Co/2 if 

(ll- 7 ) ||MI|ro+4,oo < 1 

where r = [n/2] + 1 is the Sobolev exponent. Here u = x — xq and xq is the approximate solution. 
Proof. It follows from Lemma 11.2 that 

(11.8) \\\h\\\ 2 ,oo < C 2 ( Cl , \\\u + Xo|||ro+3,OD, |||^o|||7-o+3,oo) 

In view of (11.2), the physical condition with cq replaced by cq/2 hold if T is so small that (11.5) hold 
and 

(H-9) Tc 1 C 2 (2c 1 , |||u + x ||L+3,oo, HIML+3,00) ^ c o/ 8 □ 

Recalling again that in the Nash-Moser iteration we will only consider u for which (11.7) hold. From 
now on we will therefore assume that < T < 1 is fixed and so small that (11.6) hold. 

12. Tame estimates for the inverse of the 
linearized operator in terms of the coordinate. 

We have 

Theorem 12.1. Suppose that T > is so small that the conditions in Lemma 11.3 hold. Suppose also 
that x = u + xq and J<3? are smooth in [0, T] x f2 and that 5& and u vanish to infinite order as t — > 0. 
Then there are constants K, depending on the approximate solution (x ,h ), on (cq,C\) and on r, such 
that there is a smooth solution 5x of 

(12.1) <!>'(x)5x = 5<f>, in [0,7] xU 6x \t=o = Sx \t=o = 

satisfying 



(12.2) ||<Sx| 



r 

\Sx\\ r < |||x||| r+rn+ 4_ Si00 / ||(5$|| s dr 



for r > 1 if 

(12- 3 ) ||Mlk+4,oo < 1 

Moreover 



(12.4) 



Proof. First show existence the equation where the vector field is expressed in the Lagrangian frame, 
W a = (dy a /dx i )6x i and F a = (dy a / 'dx^SQ* : 



(12.5) L W = F, in [0,T] x fl W\ f Q W ; 

and that it satisfies the estimate 



(12-6) \\W\\ r + \\W\\ r < ^^|||x||| r+ro+4 - s ,oo / \\F\\ S d T 



s=l 



for r > 1 if We have already proved this for L\W = F in Theorem 9.1, using Lemma 10.1 and Lemma 
11.3. Therefore it remains to prove the result for LqW = L\W — B^W = F, where B% is given by 
(2.63). We have ||-B 3 V7|| S < KJ2t= \\\x\\\s+r +3-k,co\\W\\ k < KJ2t=i IIMIIs+r +4-fc,oo||W / ||/ c , if s > 1. 
Applying the theorem to the equation LiW = F + B 3 W and using interpolation we get that for r > 1 

(12-7) \\W\\ r + \\W\\ r <K^\\\x\\\ r+r0+A - StOO f (\\F\\ S + \\W\\ s )d T 

8 = 1 J0 

If we put up an iteration L\W k+1 = F — BsW k , for k > and W° = then (12.6) is going to be true 
with W in the right hand side replaced by and in the left by W k+1 . It is easiest to first show 

convergence to a smooth solution and then afterwards prove the estimate (12.2). To show convergence we 
can just consider the estimate (12.6) where we include the norms of x in the constants and estimate the 
lower order norms by higher order norms. Let W k+1 = W k+1 — W k . Then L\W X = F and LiW k+1 = 
-B 3 W k , for k > 1. Hence if E k = £j =1 ||ZW J '||r + ||V^'|| r we have E k+1 < C r f* (\\F\\ r +E k ) dr. Using 

a Gronwall type of argument one therefore get uniform bounds E k < C' r J * ||-F|| r dr, for < t < T < 1. 
This proves convergence of to a smooth solution. W . Once we have a smooth solution it will satisfy 
the estimate (12.7). By a Gronwall type of argument and induction as in section 8 it follows that the 
solution also satisfies the estimate (12.6), for some other constant K. 

Finally, we want to deduce the estimate for 5x and <5<E>. The estimate (12.6) is in terms of W a = 
5xidy a /dxi and F a = 5&idy a /dx z , turning them into an estimate for 5x and 5& just produces lower 
order terms of the same form: 
(12.8) 

\\6x\\ r + \\Sx\\ r < K(\\W\\ r + \\W\\ r + \\ \x\\ \ r+2 ,oo \\ W\\ ) , \\F\\ r < K(\\6$\\ r + |||x||| r+ i )00 ||(5$||o), 

where we used the interpolation inequality in Lemma 5.7. By (12.6) ||IU||o < ||^||i < K \\F\\i dr < 
K f* || (5$ ||i dr and by interpolation 

(12-9) || |x|| | r -|_ r0 _|_4_ SjOO || |x|| | s _|_i j00 < K\\ \x\\ | r _|_ rn _|_3 j00 , 1 < s < r 

Hence (12.3) follows. □ 

If we also use Sobolev's lemma, r = [n/2] + 1, and estimate the integrals by the L°° norms and 
use interpolation and (12.2) we get 

(12.10) lll^llkoo <k{¥H r+r n ,oo + |||^|||r+2r +4,oo||^^|||o,oo)) r > 
Furthermore we want to turn it into estimates for 

(12.11) 4>(u) = $(u + x ) - $(x ) 

see (2.23). Then &'(u) = &'(u + xq). Let ip(u) denote the right inverse of <&'(u). Then since |||x||| ri00 < 
||| - u|||r,oo + |||xo|||r,oo we can again include |||xo||| r ,oo in the constants. Hence we have proven that 
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Theorem 12.2. Suppose that T > is so small that the conditions in Lemma 11.3 hold. Suppose also 
that x = u + xq and g are smooth in [0, T] x f2 and that g and u vanish to infinite order as t — > 0. Then 
there are constants K, depending on the approximate solution (xo,ho), on (co,ci) and on r, such that 
the linearized operator <J>'(it) has a right inverse ip(u) satisfying 

(12.12) III^Hffllkoo < ^(lllfflllr+rcoo + || Ml |r+2r +4,oo \\fj\\ |o,oo) , T > 

*/ 

(12-13) ||MIL+4,oo < 1 

13. Tame estimate for the second variational derivative. 

We now first want to show that the Euler map $>(x) given by (2.20)-(2.21) is C k , i.e. that <J?(x) 
depends smoothly on k parameters if x does. To be more precise, with B h = {r G R fc ; |r| < 1}, we 
want to show that $(x) — $(x ) G C°°(B k , Cqq) if x — x G C°°(B k , Cqq), where x is the approximate 
solution satisfying (2.18) and Cqq is given by (2.19). That D\ and di in (2.20) depends smoothly on 
parameters is obvious so we only need to prove that h = *(x), given by (2.21) does. Subtracting off 
the approximate solution Hq of (2.21) we get (10.2), where the right hand side is in C k (B k ,C§q) if 
x-x e C k {B k , C§§). Hence it follows from Lemma 7.6 that h - h G C k (B k , C§§). 

We must now also obtain tame estimates for the second variational derivative. If x depends smoothly 
on the parameter r then the variational derivative of x is 5x(t,y) = dx(r,t,y)/dr\ Q , we can e.g. take 
x = x(t,y) + r5x(t,y). The first variational derivative $'(x) of the Euler map is given by 

(13.1) &(x)6xi = 5<5>{x)i = d<S>(x)i/dr\ r=0 = D 2 t 5xi - d k h di5x k + h'{8x), 
where Sh = h'(Sx) = ^f'(x)Sx satisfies 

(13.2) D*{c5h)-A5h = -5Ap-d kP A5x k - 2{d i d k p)d i 5x k , 5Ah = 1d k V l drfx 1 diV k -2d k V i diSv k 

and o~p\ gfl = 0, where Sv = D t 5x. 

Now, let x depend smoothly on two parameters r and s, such that d 2 x/drds = 0, and also set 
ex = dx/ds\ _ Q , e.g. x = x(t, y) + r5x(t, y) + sex(t, y). Then the second variational derivative is given 
by 

(13.3) &'(x){6x, ex)t = e5<I>(x) i = d(d^i(x)/dr\ r=0 )/d a \ a=Q . 

We have 
Lemma 13.1. 

(13.4) $"(ex, Sx)i = d k h (d l ex l di5x k + drfx 1 diex k ) - d k h'{ex) diSx k - d k h'(5x) ^ex fc + ^ h"(ex, Sx) 

where h' = *'(x) and h" = V"(x). 

The estimates for h! = ^'(x) and h" = fy"(x) must also be obtained: 
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Lemma 13.2. Let h = \l/(x) and let 5h = h'(5x) = ^f'(x)Sx be the variational derivative. We have 

(13.5) p/i||| ri00 <tf(Px||| r+ro+l,oo + || I^H |r+2r +3,cx) || I^H |o,oo) 
and with h"(Sx,ex) = ^"(x)(5x,ex) the second variational derivative, we have 

(13.6) \\h"(Sx,ex)\\ ri00 

< K(|||(5x||| r+ 3 ro+6>00 |||ex|||o,oo + |||ex||| r+3ro+6>00 |||fe||| 0i00 + |||x|||r+3r +6,oo|||^|o,oo||MI|o,oo) 

Proof of Lemma 13.2. The proof is similar to the proof of Lemma 10.1. We have 

(13.7) D 2 t {ch) - Ah = (diV'XdjV*), Ah = k~ 1 d a {Kg ab d b h) 
It follows that 

(13.8) D 2 t {c5h) - ASh = S^V^d^)) + k' 1 d a (5(Kg ab )d b h) - div5x k' 1 d a (Kg ab d b h) , 

since 5k = KdivJx, see [LI]. Using the estimate for h in terms of x in Lemma 10.1 and Corollary 7.5 
gives if / denotes the right hand side of (13.8) 

( 13 - 9 ) lll/lllr.oo < ^(lll^ll|r+2,oo + IIMH lll fe lllo,oo) 

and hence by Corollary 7.5 

(13.10) III^IHr.oo < ^(|||^|||r+ro+l,oo + || Ml |r+2r +3,cx> || |&c|| |o,oo) 

In the proof we use the interpolation inequalities in Lemma 5.7 and the fact that || |x|| |r +4,oo < K. 

To calculate the second variation we apply e to this, where we assumed that e5x = 0. Note that 
edivJx = —(diex k )(dkSx l ). 

(13.11) D 2 t {ce5h) - Ae5h = e5((diV j )(d 3 V 1 )) 

-&\v€XK- l d a {5{ng ab )d b ti) -dw5xK- 1 d a (e(Kg ab )d b h) 
+ (divex divJx + {diex^dkdx 1 )) k' 1 d a {Kg ab d b h) - divJx k' 1 d a (e{Kg ah )d b ti) + k' 1 d a {e5{Kg ab ) d b h) 
+ K~ 1 d a (S(Kg ab )d b eh) -&\v5xK- l d a (ng ab d b €h) + n~ l d a {e{ng ab )d b 5h) - &\v ex k' 1 d a {ng ab d b 5h) 

Here 5g ab = -g ac g bd 5g cd , 5g ab = 5 ij (d a 5x l )(d b x : >) + 5 ij (d a x l )(d b 5x : >) so since ebx = 0, we have eSg ab = 
25ij(d a 5x t )(d b ex : '). The first terms on the right of (13.11) gives raise to term of the form 

(13.12) (dev)(ddv), (dev)(dSx)(dv), (ddv)(dex)(dv), (dex)(d5x)(dv)(dv) 
multiplied by powers of dy/dx or dx/dy. If /i denotes any of these terms then 

(13.13) |||/i||| r ,oo < i^(|||5x||| r+4>00 |||ex|||o,oo + |||ex||| r+4>00 |||Jx|||o,oo + II |r-+5,oo || l^^ll |o,c3o || |ea?|| |o,oo) 
The terms on the second and third row in (13.11) gives raise to terms of the form 

(13.14) (d 2 5x){dex){dh), (d5x){d 2 ex)(dh), (d5x)(dex)(d 2 h), (d5x){dex){dh)(d 2 x) 
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multiplied by powers of dy/dx or dx/dy. These can be estimated by (13.13) plus 

(13.15) i^|||/i||| r+4>00 |||(5x|||o,oo|||ea;|||o,oo < if |||x|||r+r +5,cx>|||Mlo,oo||MI|o,oo 

by Lemma 10.1. The terms on the last row in (13.11) gives raise to terms of the form 
(13.14) 

(d 2 5x)(deh), (d 2 €x)(d5h), (d5x)(d 2 eh), (dex)(d 2 5h), (dex)(d5h)(d 2 x), (d5x)(deh)(d 2 x) 
multiplied by powers of dy/dx or dx/dy. These can be estimated by 
(13.15) 

K(|||(5x||| r+ 3 )00 |||e/l|||o,oo + |||ea;|||r+3,oo|||^|||o,oo + |||<^|||r+3,oo||M||o,oo + || \etl\\ | r+3 ,oo || \Sx\\ | ,oo) 

+ #111 

x\\ |r+4,oo 

If we also use (13.5) we see that this can be estimated by 

(13.16) 

K{\\\6x\\\ 

r+3,oo || | £X 1| | ro + l,oo + || | £X \ \ |r+3,oo |||fe|||r + l,oo + |||^||| 'r+'ro+4,oo || I^H |o,oo + || |r+ro+4,oo lll fa lllo,oo) 

+ #111 •^11 |r+4,oo 

(infill 

0,oo II l^^ll |ro + l,oo + |||fe||L + l,oc||M||o,oc) 
+ #11 kill 2r n + 3,oc(|||&£|||r+3,oo|||ea ; |||o,oo + || Mil |r+3,oo || |<M| |o,oo) 

+ #111 x\\ |r+2ro + 6,oo (|||<fa|||o,oo||MH|o,oo + ||MH|o,oo|||<fa|||o,oo) 

+ ^||kl||r+4,oo||k|||2ro+3,oo(|||fe|||o,oo|||ea;|||o,oo + || Mil | ,oo || Mil | 0>O o) 

Using interpolation again this is bounded by 

(13.17) K(|||fe||| r+ro+4>00 |||ex|||o,oo + ||M|| r+ro+4>00 |||(fa|||o,oc + ||M||r+2r + 6,oo|||Mlo,oo||M|||o,oo) 

It follows that 

(13.18) \\\D 2 (ce5h) - Ae5/i||| r>co 

< K(|||5x||| r+rn+ 4 >00 |||ex|||o,oo + |||ex||| r+ro+4>00 |||fe||| 0i00 + |||x||| r +2r +6,oo||Mlllo,oo||Mlllo,oo) 

Hence by Corollary 7.5 : 

(13.19) |||e^||| ri00 

< if (||MII|r+2r +3,oo||Mlllo,oo + || Mil | r+ 2r +3,cx) || \Sx\\ | ,oo + || Ml |r+3r +5,oo || \$x\\ | ,oo || Mil | ,oo) 

Theorem 13.3. Suppose that T > is so small that the assumptions in Lemma 11.3 hold. Then there 
are constants K depending on the approximate solution (x ,h ), on (co,ci) and on r such that 

(13.20) \\\&'{u + x Q ){ex,8x)\\\ Tt00 

< if (|||<5x||| r+2 ro+4,oo||M||| ,oc + || Ml |r+2r n +4,oo || M|| |o,oo + || Ml | r+3ro+6 ,oo || \5x\\ | ,oo || Mil | 0>O o) 

if 

(13.21) |||u||| rn+ 4 >00 < 1. 
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14. The smoothing operators. 



We will work in Holder spaces since the standard proofs of the Nash-Moser theorem uses Holder 
spaces. The Holder norms for functions defined on a compact convex set B G R 1+n are given by, if 
k < a < k + 1, where k > is an integer, 

,,,,, \d a u(t, y) — d a u(s, z)\ . . 

(14.1) u a>00 = \\\u \\\ H a = sup y -T— h sup \u{t,y)\ 

{t, y ),{ s ,z)eB \(t,y)-(s,z)\»-k (tiW)eB 

and = su P(t,y)e_B l n (^2/)l- Since we use the same notation for the C k norms, |||u|||fe i00 we will 

indicate the difference by simply using letters a,b,c,d,e, f etc for the Holder norms and i,j,k,l, .. for 
the C k norms. Since a Lipschitz continuous function is differentiable almost everywhere and the norm 
of the derivative at these points is bounded by the Lipschitz constant, we conclude that for integer 
values this is the same if the L°° norm of d a u for \a\ < k, and furthermore, since all our functions 
are smooth it is the same as the supremum norm. Our tame estimates for the inverse of the linearized 
operator and the second variational derivative are only for C k norms with integer exponents. However, 
since |||u|||fc )00 < C|||u||| a)00 < C|||u|||fc+i )00 , if k < a < k + 1, see (14.2), they also hold for non integer 
values with a loss of of one more derivative. 

In [L3] we used smoothing only in the space directions but here we will use smoothing also in the 
time direction. Therefore we define the Holder space time norms as above. 

They satisfy 

(14.2) HMIkoo < C||klll&,oo, a<b 
and they also satisfy the interpolation inequality 

(14-3) IIMIIcoo^CIIMII^IHulllJ^ 

where < a < c < b, 0<A<1 and c = Xa + (1 — X)b. Furthermore, the Holder spaces are rings: 

(14.4) infill a,oo 

<C(\\\u\\\ IIMIkoo + IIMIkoollMlkoc) 

For the Nash-Moser technique, apart from tame estimates one also needs smoothing operator S$ 
that satisfy the following properties with respect to the Holder norms: 

Lemma 14.1. Let \\\u\\\ a denote the Holder norms in (14-1) with B = [0, T] x Q, where T < 1. Let 
Cqq = Cqo([0,T] x f2) be as in (14-11)- Then there is a family of smoothing operators Sg : Cqq — > Cqq, 
1 < 8 < oo such that 

(14.5) |||5Hlla > oo<C|||u||| 6 , 0O , a<b 

(14.6) |||5fl«|lkoo < Cr- b |||«||| 6i00 , a>b 

(14.7) \\\(L - S e )u\\\ a < C9 a - b \\\u\\\ b>00 , a<b 

(14.8) |||0S 2e -SeMlkoo <Ce a - b \\\u\\\ b>00 , a>0 

where the constants C only depend on the dimension and an upper bound for a and b. 

The last property, (14.8) follows from (14.6) for a > b and from (14.7) for a < b. Alternatively, it 
follows from the following stronger property 

(14.9) |||^^«|lkoo<Ce°- 6 - 1 |||«||| 6 , 00 , a>0. 
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For functions supported in the interior of a compact set K there there are smoothing operators, 
see [HI], that satisfy the above properties (14.5)-(14.9), with respect to the Holder norms. These are 
constructed as follows. Let the Fourier transform <fi £ Cq° be 1 in a neighborhood of the origin and 
set <f>e(z) = 1+n 4>(0z), and Squ = x4>e * u, where x G Co° is 1 on a neighborhood of K. However we 
have functions defined on the compact set [0, T]x£l that do not have compact support in £1. Therefore 
we need to extend these functions to have compact support in some larger set, without increasing the 
Holder norms with more than with a multiplicative constant. There is a standard extension operator 
in [S] that turns to have these properties, see Lemma 14.2 below. 

First however, we note that we will only apply the smoothing operators to functions that vanish 
to all orders as t — > 0. Hence we can extend these functions to be for t < without changing the 
Holder norm. Then we extend the functions defined for t < T to functions supported in [0,2], using 
the extension in Lemma 14.2, for y G f2 fixed, noting that Holder continuity in (t,y) follows from 
differentiability and Holder continuity in each direction using the triangle inequality and the linearity of 
the extension operator in Lemma 14.2. Then we want to extend the functions defined in Q = {y; \y\ < 1} 
to functions supported in {y; \y\ < 2}. In order to do this we first remove a region around the origin 
and introduce polar coordinates r > and uj £ S*™ -1 , which is a nonsingular change of variables away 
from the origin. Then we use the extension operator in Lemma 14.2 to, for fixed t and u extend in the 
radial direction from function defined for r < 1 to functions supported in r < 2. By the remark above, 
Holder continuity in (t, r, uj) follows from Holder continuity in each direction. 

Doing the extensions above we hence obtain an extension u of u defined in [0, T] x such that 

(14.10) |||«|||o,oo<C|||«||| 0i00 , a>0 supp(n)G{(t,y);0<t<2, \y\ < 2} 
for u in 

(14.11) C^° = C£([0,T] x T1) = {u e C°°([0,T] x H), D k t u\ t=0 = 0, k > 0} 

We note that, in fact the constant in (14.10) is independent of T. 

Once we have the extension operator we can use the smoothing operators in [H1,H2], defined for 
compactly supported functions, applied to the extension of our function. Let us call the smoothing 
operators defined in [H1,H2] Sg. These satisfy the properties (14.5)-(14.9). By (14.10) the smoothing 
operators Sgu given by the restriction of Sgu to [0, T] x Q then also satisfy the properties (14.5)-(14.9) 
if u is in (14.11). However, S$u is not in (14.11) anymore. In our estimates we will only apply the 
smoothing operators to functions that vanish to all orders as t — > and in our estimates we need also 
Squ to vanish to all orders as t — > 0. We therefore have to modify our smoothing operators so that this 
is true. Let %(i) G C°° be a function such that %(t) = 0, when t < and x(t) = 1, when t > 1, and let 
X e(t) = xiPt). Then 

(14.12) S e u = xeS e (u)\ {0T]yIi 

is in (14.11) and we claim that for functions u in (14.11), (14.5)-(14.9) hold. This follows from Lemma 
14.3 below, since the smoothing operators defined in [H1,H2] are convolution operators of the form in 
Lemma 14.3. 

Lemma 14.2. There is a linear extension operator £xt : H a ((— oo, 0]) — > H a ({— oo, +oo)), where H a 
are the Holder spaces, such that £xt(f) = f when r < 0, and 

(14.13) ||£^(/)lkoo<C1|/|| ai00 

53 



Here C is bounded when a is bounded. Furthermore, if r > —c in the support of f , where c > 0, then 
r < c in the support of Sxt(f). 

Proof. Let £xt(f)(r) = f(r), where f(r) = f(r), when r < 0, and 

/CO 
/(r-2Ar)V»i(A)dA, r>0 

where ipi is a continuous function on [1, oo), such that 

/CO pOO 
Y>i(A)dA = l, J \ k 4> 1 (\)d\ = 0, k>0, |V>i(A)| < C N (1 + X)- N , N>0 

The existence of such a function was proved in [S] where the extension operator was also introduced. 
In [S] it was proven that this operator is continuous on the Sobolev spaces but it was not proven there 
that it is continuous on the Holder spaces so we must prove this. 

First we note that if / € C k then the extension is in C k . In fact 

/CO 
fU)( r _ 2Ar)(l - 2X) j Vi(A) dX, r > 01 

From the continuity of d{f and (14.14)-(14.15) it follows that lim r ^ +0 f {j) {r) = f ij) (0), that / is in 
C k , and that ||/|U,oo < Cfc||/||fc,co, if k is an integer. 

Suppose now that k < a < k + 1 where k is an integer. We must now estimate 

(MAT, sup -/ { y 

r,p \r-p\ a ~ k 

by Ck\\d k f\\ a _k,oc- If r < and p < there is nothing to prove. Also ifr<0<porp<0<r, then 
\f — p| > \p\ and |r — p| > \r\ so in this case, we can reduce it to two estimates with either r = or 
p = 0. Also it is symmetric in r and p so it only remains to prove the assertion when r > p > 0. It 
follows from the Holder continuity of f^ and the last estimate in (14.15) that for r, p > 0, 



(14.18) J (f( k )(r-2Xr)-f( k \p-2Xp))(l-2X) k M^dX < C k \\fW \\ a _ k>OQ \r - p\ a ~ k 
which proves the lemma. □ 

Lemma 14.3. Let \(t) £ C°° be a function such that %(i) = 0, when t < and x(t) = 1? when 
t > 1, and let Xe(t) = x(0i)- Let the Fourier transform <f> G be 1 in a neighborhood of the origin 
and let Xi ^ Co° be 1 on a neighborhood of {(t,y); < t < 2, \y\ < 2}. Set Sg = Xi4>e * u > where 
M{t,y)) = (f>(0{t,y))/e 1+n . Then 

(14.19) |||(1 - X e)S e u\\\ a ^ < C|^r- & |||n||| fe>00 , 
if u is smooth and vanishes for t < and for t>2. 

Proof. First we note that by interpolation it suffices to prove the estimate for a = k an integer. Since 
u vanishes to infinite order as t — > we have iik<b<k + l 



^ / /".6-*„,..„, 



(14.20) \u(t,y)\= / (dtu)(8,y) \, ds 



< I s ~' c ||k||| b;OO v " ds < C b t b \\\u\\\ b 



(k-iy. -Jo B,ou (k-iy. 

Since 4> is fast decaying we have if \a\ = k, 

ds dy 

/ (1 + \6t - s\ + \6x-y\) N 
for any N. Here the integral is uniformly bounded when 9t < C so the lemma follows. □ 
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(14.21) \D a {(l-xe)xi<Pe*u)\<C N e k - b \\\u\\\ b , c 



15. The Nash Moser Iteration. 



At this point, given the results stated in sections 11-14, the problem is now reduced to a com- 
pletely standard application of the Nash-Moser technique. One can just follow the steps of the proof of 
[AG,H1,H2,K1] replacing their norms with our norms. The main difference is that we have a boundary, 
but we have constructed smoothing operators that satisfy the required properties for the case with a 
boundary. The proof of the Nash-Moser that we outline below is similar to the one in [L3] . The only 
difference is that now we also smooth in time. We will follow the formulation from [AG] which however 
is similar to [H1,H2]. The theorem in [AG] is stated in terms of Holder norms, with a slightly different 
definition of the Holder norms for integer values. However, the only properties that are used of the 
norms are the smoothing properties and the interpolation property in section 14, which we proved with 
the usual definition, i.e. the one used in [HI]. 

Let us also change notation and call l>(n) from last section $(tt). For k < a < k + 1, where k > 
is an integer, let 

/ 1K1 x \d a u(t,y) - d a u(s,z)\ 

(15.1) ||M||o,oo = |||«|||fr» = sup 2^ [7T1T7 \\a-k + sup J u (*>2/)l 

(t,y),(s,z)€[0,T]xn {al = k WiV) y S > Z >\ (t,y)€[0,T]xn 

and IIMHho = sup( t y ^ € [ t]xTl V)\- The estimates we proved for the inverse of the linearized operator 
and the second derivative of the operator where in terms of L°° norms, i.e. Holder norms for integer 
values. However, since |||u|||fc )00 < |||«||| a ,oo < IIMIIfc+ii if A; < a < fc + 1, it follows that they also hold 
for non integer values with loss of an additional derivative. 

(Hi): $», is twice differentiable and satisfies 

(15.2) |||$"(u)(t;i,U2)|||„,oo 

< C a (j| \vi 1| \a+n,oo || \V2 || U',oo + || \vi \\ \fi',oo \\ \V2 1| \a+fi,oo + || Ml |a+At,oo || \vi 1| \fi',oo \\ \V2 || U',oo) , 

where (i = 3r + 8, for u, v±,V2 S Cgg, if 

(15.3) |||u|||/i',oo < 1, // = r + 4 

where K is the constant in Lemma 12.1 

{Ti.2)'- If u € Cqq satisfies (15.3) then there is a linear map ip(u) from Cqq to Cqq such that 
&(u)ip(u) = Id. It satisfies 

(15.4) |||^u) 5 ||| 

a,oo <C a III fflllo+A.oo + IHS'lllo.oo ||MIU+A,oo)) 

where A = 2r + 6. 

Theorem 15.1. Suppose that $ satisfies (Hi), (H2) and $(0) = 0. Suppose that \i > [i! and let 
a>\ + [i + [i',a^N. Then 

i) There is neighborhood Ws = {/ € Cqq; |||/||| a +A,oo < $ 2 }, S > 0, such that, for f G W5, the 
equation 

(15.5) = / 
has a solution u = u(f) £ Cgg. Furthermore, 

(15-6) |||«(/)|||o,oo<C|||/||| Q+Ai00 , a<a 
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In the proof, we construct a sequence Uj £ Cqq converging to u, that satisfy |||%|||^' < 1 and 
IH'S'i^illU' < 1) f° r an J) where 5$ is the smoothing operator in (15.7). The estimates (15.2) and (15.4) 
will only be used for convex combinations of these and hence within the domain (15.3) for which these 
estimates hold. 

Following [H1,H2,AG,K1,K2] we set 

(15.7) u i+ i = Ui + Sv,i, Sui = ip(SiUi)gi, u = 0, Si = S 0i , 6 i = 6 2' 1 , 6 > 1 
and gi are to be defined so that Uj formally converges to a solution. We have 

(15.8) - $(ui) = &{ui){u i+1 - Ui ) + e'! = '(ui)ip(SiUi)gi + e'( 

= (&{ Ui ) - ^(SiUi^SiU^gi + e'( = ^ + e'( + 9i 



where 






(15.9) 




- $'(SiUi))5ui 


(15.10) 






(15.11) 


i i a 

C{ 6^ "T 6j 




Therefore 







(15.12) - = e, + 5i 

and adding, we get 



(15.13) *(t*i) = 2* + ^ + e * + ( J " 5 0^i> ^ = E e ^ 

3=0 j = l 

To ensure that <&(iii) — ► f we must have 

i 

(15.14) + = 
Thus 



(15.15) #o = So/, g i = (S i -S i - 1 )(f-E i - 1 )-S i e i - 1 

and 



(15.16) = + et + (I - Si)Ei 

Given uo, u\, these determine Suo, 5u\, Su^ which by (15.9)-(15.10) determine ei,...,ej_i, which 
by (15.15) determine gi. The new term Ui + \ is the determined by (15.7). 
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Lemma 15.2. Assume that |||Mi||| M ',oo < 1, |||^i+i|lk,oo < 1, |||5i«i||| At ' i0 o < 1- TTten 

(15.17) |||e-||| r)00 < C r (|||(7 - Si)uj || | r+At)OC III^H l^oo + |||(I - S , i )u i ||| M ' j00 |||(5'Ui||| r . +M ) 

+ C r || | SiUi || Ir+^oo || I {I ~ SijUi || |^',oo || \SUi || |/i',oo) 

anc? 

(15.18) ll|eJ , |||r,ao<C r (|||5u i ||| III^IH /i' ,oo + llklll lll^lll 
Proof. The proof of (15.17) makes use of 

(15.19) ($'(uO - &(SiUi))6ui = [ &'{S l u l + s{I-S l )u l ){u l -S l u l ,5u l )ds 

Jo 

together with (15.2). The proof of (15.18) makes use of 

(15.20) <f>(u i+1 )-$(ui)-&(ui)6ui= [ (1- s)<f>"(ui + s6ui)(6ui,5ui)ds 

Jo 

together with (15.2). □ 

Let a > a and a — fj, > 2a — \i — //. Throughout the proof C a will stand for constants that depend 
on a but are independent of i and n in (15.21). 

Our inductive assumption (H n ) is, 

(15.21) |||«y«i|||o,oo < SO?'"' 0<a<a, i<n 

If n = then 5uq = ip(0)Sof, and if a < a, we have |||<5uo||| ,oo < Cs|||/HU+a,oo < Ca^ 2 , so it follows 
that (15.21) hold for n = if we choose 5 so small that Ca<5 < 6>o _a . We must now prove that (H n ) 
implies (H n+1 ) if C' & 6 < 1, where C~ is some constant that only depends on a but is independent of n. 

Lemma 15.3. // (15.21) hold then with a constant C a independent of i < n 

i 

(15.22) ^|||<5u J -||| Oi00 <C o <5(miii(i,l/|a-a|) + l)(^- a + l), 0<a<a 

3=0 

Proof. Using (15.21) we get E}=o lll^jllkoo < C a 5 £} =0 2j{a ~ a) and noting that Y! j= o 2 ~ sj < 
C(min (1 + l/s, i) + 1), if s > 0, (15.22) follows. □ 

It follows from (15.22): 

Lemma 15.4. If (H n ), i.e. (15.21), hold, and a > a, then for i < n + 1 we have 

(15.23) HNIkoc <C a 5(mm{i,l/\a-a\) + l)(er a + 1), 0<a<a 

(15.24) \\\S tUl \\\ a>00 <C a 5(mm(i,l/\a-a\) + l)(e^- a + l), a>0 

(15.25) |||(/-5 i )«i|||a,oo<C Wr a , 0<a<a 
where the constants are independent of n. 

Proof. The proof of (15.23) is just summing up the series Uj+i = J2)=o ^ u i-> usrn § Lemma 15.3. (15.24) 
follows from (15.23) using (14.5) for a < a and (14.6), with b = a for a > a. (15.25) follows from (14.7) 
with b = a and (15.23) with a = a. □ 

Since we have assumed that a > fi', we note that in particular, it follows that 

(15.26) |||^j|||/i',oo < 1 and |||<Si'Mi||| M ',oo < 1, for i<n + \ if C^S < 1. 

As a consequence of Lemma 15.4 and Lemma 15.2 we get 
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Lemma 15.5. // (H n ) is satisfied and a > y! , then for i < n, 

(15.27) ||K||U, 0O <C' a <5X~ (2a ~ M ~' i ' ) , 0<a<a-y 

(15.28) |||<||U,oo < CJ 2 er i2a ~^'\ < a < a - y, 

where the constants are independent of n. 

As a consequence of Lemma 15.5 and (14.8) we get 
Lemma 15.6. // (H n ) is satisfied, then for i < n + 1, 

(15.29) HI -Vi-iHkoo < C a 5 2 er (2a ~^'\ a > 

(15.30) |||(5 i -5 i _i)/||| 0>00 <C a (9r /3 |ll/llkoo, «>0 

(15.31) |||(/ - ^)/||U,oo < C^r^lll/llkoo, < a < (3 

Furthermore, if a — y > (2a — fi — y'): 

(15.32) 111(5, - 5 i _i)£; i _ 1 ||| , 00 < C a S 2 er i2a ~^'\ a > 

(15.33) \\\{I-S i )E l \\\ a ^<C a 5 2 e a l - {2a ~ ,1 ~ tl '\ 0<a<a-y 
Here the constants C a are independent of n. 

Proof. (15.29) follows from (15.27). For a < a — y we use (14.5) with b = a and for a > a — y, we use 
(14.6) with b = a-fi. (15.30) follows from (14.8) and (15.31) follows from (14.7). Now, Ei = Ej=o e j 
so by Lemma 15.5 |||.Ej|||,5_ M!00 < C a 5 2 Y?j=o 0^~^~ ( - 2a ~ fl ~^ ' < (j' a g2ga-n-{2a-^-n ) ^ s ^ nce we assumed 
that the exponent is positive. (15.32) follows from this and (14.8) with b = a — \x and similarly (15.33) 
follows from (14.7) with b = a — \i. □ 

It follows that: 

Lemma 15.7. // (H n ) is satisfied, a — \x > (2a — fi — y!), and a > y! then for i < n + 1, 

(15.34) |||ft|||a,oo < CjHr^-^ +C a 0r /3 |ll/lll/3,oo, a > 0. 

Using this lemma and (15.4) we get 
Lemma 15.8. // (H n ) holds, a — \i > (2a — y — pi!), a > fi' , a > X then, for i < n + 1, we have 

(15.35) Hl^llkco < c^e^-^-^ + C a |||/||| /3 , 00 ^ +A - /3 , a > 0. 
Proof. Using Lemma 15.7, (15.24), and (15.4) we get 

(15.36) IllJUilH^oo < C a (|||ft||| a+ A !00 + |||5||U,oo|||'S , i^|||a+A,oo) 

< Ca (5 2 9^ x - {2a -^ + |||/||| /3 , 0O C +A " /3 ) 

+ C a (S 2 9t [2a -^' ] + Hl/llkoo ^^minCi, l/\a - a - A|) + 1M +X ~ a + 1) 
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Using that a > d we get (15.35). □ 

If, we now pick (3 = a + X, and use the assumptions that X + a < 2a — fx — //, and |||/|||ch-a,oo < $ 2 : > 
we get that for i < n + 1, 

(15.37) lll^llkoc < C a 5 2 e?- a , a > 0, 
If we pick 5 > so small that 

(15.38) C & 5 < 1, 

the assumption (i? n+ i) is proven. 

The convergence of the is an immediate consequence of Lemma 15.3: 

oo 

(15.39) ^ || \u i+ i - Ui\\\g t00 < C a S, a < a 

i=0 

It follows from Lemma 15.6 that 

(15.40) WMm) - /nUoo < c a 5 2 er a ~ x 

which tends to 0, as i — > oo, if a < a + A. In particular it follows from (15.39) that IH^IIU'.oo — 1> ^ 
5 > is chosen small enough. 

It remains to prove u G Cgg. Note that in Lemma 15.8 we proved a better estimate than (H n ). In 
fact if we let 7 = 2(a — //) — (a + A) > and a' = a + 7, then |||/||| q '+a,oo < C implies that 

(15.41) |||<5«i||koo<C o er a, ) «>0 

Using this new estimate, in place of (H n ), we can go back to Lemma 15.4-Lemma 15.8 and replace a 
by a' . Then it follows from Lemma 15.8 that 

(15.42) lll^llkoc < Ca e a t +x - 2{a '-^ + cx^wmUoo 

and if we now pick 7' = 2(a'— //) — (A— a') = 27 and a" = a'+^f' = a+27, and use that |||/||| a '+ 7 ',oo < C 
we see that 

(15.43) |||<5«i|||o,oo<C o e?- a ", «>0 

Since the gain 7 > is constant, repeating this process yields that (15.41) holds for any a' and hence 
that (15.39)-(15.40) hold for any a > 0, (without S ). It follows that uj is a Cauchy sequence in 
C fc ([0,T] x fi), for any k, and since Uj G it follows that u,- — > u G Cqo and — > / G Cgg, and 
since $ is continuous it follows that $(tt) = /. (15.6) follows from (15.37) with <5 2 = |||/||| a +A,oo- This 
concludes the proof of Theorem 15.1. 
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16. Existence of initial data satisfying the compatibility conditions. 



In this section we show that there are initial data satisfying the compatibility conditions. We do 
not attempt to find the most general class of initial data that do so, the purpose is simply to show that 
our local existence theorem is not about the empty set. The set of initial data we construct can then 
easily be extended to a much larger set using essentially the same proof. 

Let us therefore start by making some simplifying assumptions. First we assume that e(h) = h. We 
now want to find a formal power series solution in t of the system 

(16.1) D^Xi = -dih, D 2 t h -Ah = (diV'XdjV*), Vi = D tXi 
with initial data 

(16.2) x\ t=0 = f , D t x\ t=Q = v , h\ t=Q = h , D t h\ t=0 = hi = -dwV 

Lemma 16.1. Let hi = D\h\ t=Q , for I > and let Vq = v\ Q . If (16.1) and (16.2) holds then 

(16.3) hi +2 = Ahi + F l+2 (h- 1 ,...,hi- 1 ), l>0 

where F 2 = (c^Vq )(e?j-Vo) and in general Fi +2 = Fi +2 (h-i, hi-\) is a sum of the form 

(16.4) F l+2 = C} x \± n h%---h%, h? = d:h h l>0, h a _ 1 = d«'V^, a=(a',i k ) 
with 

(16.5) |cki| + h + ... + \a n \ + l n = I + 2, n > 2, -1 < l { < I - 1, 1 < + l { < I + 1 

Proof. If we us that [D t ,di] = — (diV k )d k we obtain from differentiating the first equation in (16.1) 

(16.6) diD l t h = -D l t +1 Vi + a^-^^VfcJ • • • (dV^V^ 

where = D^V . Here the sum is over k\ + ... + k n = I + 2 — n, n > 2, k n > 1, and terms in the sum 
consists of contractions over n — 1 pairs of indices. From differentiating divV = diV 1 , 

(16.7) D l t +1 divV = divD l t +1 V + d lkl - k "(dV kl ) ■ ■ ■ dV kn 

where the sums are over k\ + ... + k n = I + 2 — n, n > 2 and terms in the sum consists of contractions 
over n pairs of indices. 

It follows from (16.1) that D t {D t h + d\vV) = 0. Hence 

(16.8) AD l t h - D\ +2 h = e lk ^- k -{dV kl ) ■ ■ ■ 8V kn + d lk ^(dV kl ) • • • (flV^J^,, 
where k n > 1 in the last sum. Finally we note that we can turn V k+ i, for k > 1 into h k : 

(16.9) V k+1 = -dh k + a kk ^- k "(dV kl ) ■ ■ ■ {dV kn _,)V kn 

where k n > 1 and V\ = —dh. This proves the general form (16.4) and it remains to prove the range 
of the indices in (16.5). To prove this we note that the terms in (16.8) are contractions over n pairs of 
indices and this is still true for the terms we obtain by using (16.9) to replace the factors of V by factors 
of h. This proves that |cki | + ... + |q„| = 2n. On the other hand when we replace factors of V k+ \ by h k the 
number of time derivatives go down by one for each factor, so we conclude that l\ + ... + l n = I + 2 — 2n. 
This proves (16.5) apart from the last statement. That \a.i\ > 1 is clear and if U = — 1 then \cti\ > 2, in 
view of (16.4), so in general \ati\ + Zj > 1, and since n > 2, (16.5) follows. 

We will now obtain a formal power series solution in the distance to the boundary of the system 
for hi in Lemma 16.1. In order to do this, we will first choose simpler initial data for (16.2), 

(16.10) f = y, v = d<P, k G = -A<P 

where 4> is to be determined. Let h-\ = —<j), then (16.3) hold also for / = —1 with F\ = 0. 
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Lemma 16.2. Suppose that g a b = 5 a b- Suppose also that ho,i and h\ t i are smooth for I > —1, and let 
Fi be as in Lemma 16.1, and F\ = 0. Then the system 

(16.11) A/i, = /i i+ 2 + F i+ 2(/i-i,...,/ii-i), l>-l 
with boundary conditions 

(16.12) h l \ dn = h ,i, V N h l \ m = h 1 , l 1>-1 
has a formal power series solution in the distance to the boundary: 

(16.13) ^,(r,w)~X)ViH(l-rr/n! 

Let x be smooth such that x(d) = 1; when \d\ < 1, x(d) = 0, when \d\ > 2 and x > 0. Then there are 
£i n > such that 

(16.14) h t (r, lo) = J2 X((l " OMnK, - r)"/n! 

are smooth functions and such that (16.11) hold to infinite order at the boundary. 
Proof. We have 

(16.15) A = + tr Vzv + A 

where 9 is the second fundamental form of the boundary and A is the tangential Laplacian on the 
boundary. In the case g a b = 5 a b, V/v = d r and ti9 = (n — l)/r, where r is the radial derivative. 
Furthermore A = r~ 2 A UJ , where A w is the angular Laplacian on S*™ -1 . Hence we have the system 

(16.16) d%=-^^h l + h l+2 -^d r h l +F l+2 {h- 1 ,...,h l -i), i>-i 
We want this system to be satisfied to all orders at the boundary, so if 

(16.17) fcfc,J = #>u|en> F m:l = d^Fih.u ...,h l . 3 )\ m 
we want 

(16.18) h k j= c ijkm^ljh m ,l+2j + ^2 dij jfc m A^F m)J+ 2j 

2i+2j+m<fc,m<fc-l 2i+2j+m<fc 

to hold for all /c > 2 and / > — 1, where hoj and h±j are the given boundary conditions. 

We now want to use induction. Note that the first term in the right of (16.18) contains hf-'j/ for 
k' + V < k + 1, and k' < k, and the second term contains hk>,i> for k! + V < k + 1, by the last inequality 
in (16.5). Assume that we found 

(16.19) h kJ , for k + l<N, I > -1, fc > 

such that (16.18) hold for k + I < N and k > 2. Note that if iV = then k < 1 so there is nothing to 
prove. We know want to find h k , i for k + I = N + 1 such that (16.18) hold also for k + 1 < N + 1. This 
is again proven by induction. Assume that in addition to (16.19) we found 

(16.20) h kth for k + l = N+l, and < k < M, I > -1 
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such that (16.18) hold for k > 2. Note that for M < 1 there is nothing to prove. Since m in the first 
sum on the right of (16.18) is less than k in the left it follows that we can find fofc+i, ;_i such that (16.18) 
hold. 

In order to prove (16.14) we note that h m j are smooth functions on S n ~ x . Hence we can use the 
usual trick of choosing e m i so small that (||/i m ,j|| m +i + l)e m z < 1/2, in which case the sum converges 
in H m for any m and r. □ 

Now, we want to find a formal power series solution in t of the system (16.1) with initial data in 
(16.2) of the form 

(16.21) f (y) = f (y) + y, v = v -dh- 1 , h = h + h h x = - divV + A h-i 

where ho, h_\ are given by Lemma 16.2 and fo, Vq and ho vanish to infinite order at the boundary. Let 
hi, for I > 2 be defined by (16.3). Then it inductively follows that 

(16.22) h^hi + ht, 

where hi are as in Lemma 16.2 and hi vanish to infinite order at the boundary. Therefore if we choose 
boundary data in (16.12) such that ho, i = for I > 0, hi,o < Co < 0. Then it follows that we can choose 
ho so that ho > in and Viv/io| an < — Co < 0. Moreover it follows that 

(16.23) ^L = °> ^° 

and hence the compatibility conditions are satisfied to all orders. 

We can now construct smooth functions in [0, T] x Q, satisfying the initial conditions (16.2) and 
the equations (16.2) to infinite order as t — ► 0: 

Lemma 16.3. Suppose that initial data f , vo and h and hi = — divVb satisfy the compatibility 
conditions for all orders, i.e. if hi for I > 2 are defined by (16.3) then 

(16.24) ^L = ' f° r l >°- 



Then there are smooth functions (x,h) in [0, T] x Q, such that (16.2) hold, (16.1) is satisfied to infinite 

order as t — > and h\ „„ = 0. 

I all 

Proof. Let x t> e as i n Lemma 16.2 and set 

oo 

(16.25) h{t,y)=Y J X{t/ei)hi{y)t l /l\ 



1=0 



where Ei > are chosen so that (||/ij|| m + l)s m < 1/2. Then it follows that the sum converges in H m 
for any m so h is smooth and satisfies h\ dQ = and D l t h\ t=Q = hi. Furthermore let x(t,y) be defined 
by 

(16.26) D 2 t Xi = -dih, Xi \ t=0 = f , D t x\ t=Q = v . □ 



62 



17. The general case, when the enthalpy is a 
strictly increasing function of the density. 



We will now outline how to generalize the existence result obtained for e{h) = ch to the case when 
e(h) is a smooth strictly increasing function satisfying 1/c^ < e'(h) < c\. First we will show that the 
functional h = ^(x), i.e. the solution of (2.21), exist for x in a bounded set, |||it||| T . + 7 . ljOO < 1, and for 
T sufficiently small. Since in the Nash-Moser iteration we need a bound for |||/i|||3,oo, we will first show 
that we can obtain such a bound as well as a bound for ||/i|| ro+ 3 independent of T. Local existence for 
the nonlinear wave equation follows from a standard argument using essentially the same estimates, so 
it is just a question of showing that we have a priori bounds up to some time T > that only depends 
on the approximate solution {xq,Hq) and is independent of x = u + xo as long as || | r0 _|_ rijOO < 1. 
Once we have a bound for || |fo|| |4 i00 , the bound for higher derivatives follows from this. 

The equation we study in this section is 

(17.1) D t (e'(h)D t h)-Ah = f, h\ t=Q = 0, f = (0,-^) 
where 

(17.2) e' + l/e' < c[, ^ \g ab \ + \g ab \ < c[, \dx/dy\ 2 + \dy/dx\ 2 < c[ 

a,b 

for some constant < c[ < oo. Let K' w denote a continuous function of d x that also depends on the 
order of differentiation r but is independent of a lower bound for T. 

We prove the following Theorem: 

Theorem 17.1. Let r = [n/2] + 1 be the Sobolev exponent and let k > 1. There are continuous 
function Ck and such that ifT>0 is so small that 

( 17 - 3 ) TC7i(|||x |||ro+2+l,oo, |||M|r +2+l,oo) < h 

and u = x — Xo is small that 

( 17 - 4 ) |||w|||r +2+fc,oo < 1 

then (17.1) has a smooth solution for < t < T satisfying 

( 17 - 5 ) III^IHfc,oo < -Dfc(||ko|||r +2+fc,oo, |||^o|||r +2+fc,oo)- 

Furthermore, there is a continuous function K' 2 of 

(17.6) IIMIkoo + |||/»|||2,oo + WMkoo + V^ + C'l 

depending also on r such that 

( 17 - 7 ) III^IHr,oo < ^2(lll a; lllr-+r-o+2,oo + || |^0 || |r+r +2,oo) 



Proof. Since by Sobolev's Lemma 

(17.8) IINkoo^Cdl/iHi^ + ii/Mi,,,,,,^ 
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where C is independent of T. It follows from Lemma 17.2 that 

(17.9) dE r /dt < M r (E r+1 + l) r + 2 , r > r 

where M r = ^(olll^lllr+^oo" Integrating this inequality gives -d(E r+1 + l)-( r+1 )/dt < M r (r + 1) and 
hence (E r+1 (t) + l)-( r+1 ) > (E r +i(0) + l)"( r+1 ) - M r (r + l)t. If M r (r+l)t < (^ r+ i(0) + l)-( r+1 )/2 it 
follows that E r+1 (t) +1 < 2(£ r+1 (0)+l). Since £ r+1 (0) < ^ EI=o [[H]] r _ s ^(||Mi, s + IIM(Wi) 
this proves the first part of the theorem for k = 1 and for k > 2 it follows from also using Theorem 
17.3. 

It follows from Lemma 17.2 and interpolation in space time that 

dE r 

(17.10) < K' 2 {E r+1 + |||a:||| r+ 3 i00 ), if E r+1 = ^ |||x||| r+2 - s -B s +i 

Multiplying by the integrating factor e^ 2 * and integrating, we get 

(17.11) E r+1 (t) < K 2 (E r+1 (0) + ||| 
By Theorem 7.3 ||/i|| r +i < K' 2 E r+ i and hence 

r 

(17-12) IIMt»-)||r+i <^(^|||x||| r+2 _ s ||/i (0,-)|| s+1 + HMHr+3) 

Using Sobolev's lemma and interpolation the estimate for |||/i||| r ,oo follows. □ 

Lemma 17.2. Suppose that g ab and e' = e'(h) are smooth and satisfy (7.4)- For s > let 

(17.13) E s+1 (t) = (lje'(p*+ 1 h) 2 + g ab {b s t H a ){t) s t H b )Kdy) 1/2 , H a = g ab d b h, 
E Q (t) = J n h 2 ndy and E r+1 = ££+J E s- Then forr>0 

r 

(17.14) \\h\\ hr + \\h\\ , r+1 <C^ [[[doo]]] r _ Sj00 E 1+s 

s=0 
r 

(17-15) E 1+r < [[[d*]]] r - s ,J\\h\h,s + \\h\\o,.+i) 

s=0 

and 

~ r +l min(r+l— k,r) 

("■is) ^ <^oEwU( E im]] r+1 - k - s , 00 E s+1 + [[[0 X ]]] r+2 _ k ^) 

k=0 s=0 

Proof. The proof is similar to that of Lemma 7.2. We have 

(17.17) -£± = y^(e'(^ 1 M(A S+2 M+5a 6 (A S ^ a )(A S+1 ^ b ))«^ 

+ l -f ((Ae')(A s+1 /i) 2 + (A5a 6 )(A S ^ a )(A^ 6 ))«^- 
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Here the terms on the second row are bounded by i£"{ (||/t||i iOO + ||fl'||i,oo)-£'s+i- The ^ rst row °f the 
right hand side is up to lower order 

(17.18) jj{D s t +1 h){b s t D t {e'D t h)) + {D s t H a ){d a D s t +l h)) K dy 

= J (D s t +1 h) (b s t D t (e'D t h) - K - 1 d a { K D°H a )') ndy = J {D s t +1 h)(D s t f) K dy, 

where we have integrated by parts using that D^ +1 h\ m = 0, that D t (e'D t ti) — n~ 1 d a (nH a ) = / and 
that Dl{n~ 1 d a ( y KH a ^ = n~ 1 d a {nDlH a ^. In fact, using Lemma 6.3 we get, since H_ a = d a h, 

(17.19) D s t +1 H a - g ab d a D s t +l h = - £ ( * + 1 )g ab (D 8 + 1 - i g bc )D\H c 

i=o V 1 J 

and, since D t D%e(h) = k' 1 D s t (kD 2 e(h))) , 

(17.20) D?D?e(h) - e'(h)D s t +2 h = £ ( S ) ^{D^ k)^ e{h)) + D 2 +*e{h) - e'(h)D?+ 2 h. 

t=0 w 

Here the right hand side is bounded by a constant times 

s-1 

(17.21) ^K-^r^l \e W (h)\\D?h\---\D' t »h\+ \e {k) {h)\\D^h\...\Dl"h\ 

i=0 Sl + ... + s k =i+2, Sl + ... + s k = s+2, 

S;>l,fc>l S i >l,fc>2 

The L 2 of this can be estimated by Theorem 7.3: 

s+1 min(s+2 — fc,s+l) 

(17.22) K' lfi J2 \\h\\U( J2 [[[dxm^-^JWHo, + ||A|| M -i) + [[[dx]]] a+2 _ kt0O ) 

k=0 i=l 

By Lemma 6.3 D s t {g ab d b h) = g ab d b D s t h + Y.lZl { s ^g ab {b s t ~ l 9bc)b l t {g cb d b h) so it follows that \\dD s t h\\ < 
[[[ dx ]]] r ^+* Since also llAVII < K w [[[dx]]] s+2 the lemma follows. □ 

Theorem 17.3. Suppose that D 2 e(h)-Ah = f, where f = (^^K^Y*). Suppose also that \e (k) (h)\ < 
Ck- Then we have 

(17.23) || [h] j r || < K[ fi Y [[[^]]] r _ fc - i> co(ll fc ll*-i.i + + [[ML-*,*,) 

k=j—i i=i 

and 

(17.24) || [h]i || < K[ fi Y WHUiX [[MU-.ooOHIo, + \\Hi,i-i) + [[[^]]],- fc ,oo) 

k=j — l i=l 

w/iere 

(17.25) = Y \hW--\h\ rk , [[[0x]]] riOo = lll^llk,oo---|||^||| rfci00 

ri + ...+r k =r, r-i + ...+r k =r, 

ri>l,k>j Ti>l 

Proof. The first inequality follows from Lemma 17.4 below and interpolation and the second follows 
from the first and Lemma 17.6 below. □ 
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Lemma 17.4. Suppose that D^e(h) - Ah = f, where f = (d i V J )(djV l ). Suppose also that \e^ k \h)\ < 
Cfc. With notation as in Definition 5.1 we have 

("■26) [h] r ,s<K' w £ 

i+j<r+s 

where [dx\ l is as in Definition 5.1, 

(17.27) [h] rt8 = \ h \rus.---\h\r k ,s k , [h] r ^ m) = E \ h \r^---\h\ 

ri + ...+r k =r, ri + ...+r k =r, 

Si + ... + S k = S Si + ... + S k = S 

Si+ri>l Si+ri>l, Si<m, 

and [h] Qfi = 1, [h] {01) = 1. 
Proof. If /i = 77(e) is the inverse of e(/i), then 

(17.28) Nr,.<C E 1^)11 

ri + ...+T- fc =r, si + ... + s fe = s, r; + s;>l 



e lri,si " " " \ e \r k ,s k - 



Since, 

(17.29) |e| r , s < |A/i| r>s _ 2 + |/| r>s _ 2 , s>2 
where 

(17-30) |A/i| r>s _ 2 < K 10 E H^.^ki, s>2 

l<i<r+2, j<s-2 

and 

(17.31) \f\ r ,s-2<K 10 [dx] r+s , s>2 

we obtain 

(17-32) |e| r , s <^ 10 E [fe]^.,.^^,, ^>2 

i<r+2, j<s-2 

It follows that 

(17-33) W r , ( ., m) <^o E M^-jl^M- 

and the lemma follows by induction. □ 

Lemma 17.5. For r > 1 we have with notation as in Definition 5.3, 

(17.33) ||% < KioE lKoo(E [[[^UU-i.oo (II^IIO, + ll^lll,i-l) + [[[^]]] r - h oo) 

3=0 i=l 
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Proof. We have 

(17.34) \Ah\ ryS < \e(h)\ r , a+ 2 + \f\r,s 

where 

r+s 

(17.35) \e(h)\ r!S <e'(h)\h\ r , s + Cj2\e [k) (h)\[h] r k s , where [h] r k s = E \h\ ri , Sl ■ ■ ■ \h\ rk , Sk 

fc = 2 ri -|_..._|_ r{ . =T) 

S! + ... + Sfc =S 

r; + s;>l 

and by the previous lemma and interpolation in space only 

(17.36) \\[h] r k J < K[ [[N]]H-.-i-i,ooll fc lli,co(ll^lli,o + l|fc||i-i,i) 

i+j<r+s, j>k — l,i>l 

Using Theorem 6.1 we get for r > 2 

(17.37) ||^||. +2 , s < ^T 10 ^ [[[^]]]. +s _,_ J , 00 II^IU,i+2 + ^io^ [[[^Jll^+^^ll^lli^oo 

i<r, j<s j=0 

+ *io[[H]] r+a+2>00 + *io E [[[^]]], +s+2 - l - J ,ooll /i llloo(ll^ll,0 + 

i+j<r+s+2, i,j>l 

Hence using induction, we get 

(17.38) \\h\\ r <K[ £ [[[«x]]] f ._ if00 (l|fc||o,i + ||fc||i,i-i) 



l<i<r 

^ rrr 



+ ^io[[[<k]]U + *'io E [[[ & 

Since ||/i||i < ||/i||o,i + ||^||i,0 the lemma now follows from induction. □ 
Alternatively, one can use interpolation interpolation in space-time. 
Lemma 17.6. We have 

(17.39) 111^111^-0 < CTmi^iiiDrhiii 1 - 1 ^ 

\ k £2\\\D?h\\\» 



I T — % — J , OO 1 



ll,ooll"ll' 



(17.40) HI \D\^h\ ■ ■ ■ \D l t *h\ \\\ L , < CTlll/illH- 1 !!!^' 



where m = l\ + ... + Ik- Furthermore, if h = h + h where h vanish to infinity order as t — > then, with 
a constant independent of T but depending on h we have 

(17.41) (I) \D\^h\ ■ ■ ■ \D\*h\ \\\ L , < C m {\\\h\\\^ |||/ l o|||oc, m )(E HI^HU 2 + !) 

j=o 

Proof. The interpolation inequalities are standard and it is also standard that the constant is in- 
dependent of T if h vanish to infinite order as t — > 0. Hence in the product we can estimate 
\\\D l t l h\\\ L 2/( m -i t) < IWD^hWltf/im-ii) + IIID^/ioIIIl 00 - □ 

Acknowledgments. I would like to thank Demetrios Christodoulou, Richard Hamilton and Kate 
Okikiolu for many long and helpful discussions. 

67 



References 



[AG]. S. Alinhac and P. Gerard, Operateurs pseudo-differentiels et theorem de Nash-Msoer, Inter Editions and CNRS, 
1991. 

[BG]. M.S. Baouendi and C. Gouaouic, Remarks on the abstract form of nonlinear Cauchy-Kovalevsky theorems, 
Comm. Part. Diff. Eq. 2 (1977), 1151-1162. 

[CI]. D. Christodoulou, Self- Gravitating Relativistic Fluids: A Two-Phase Model, Arch. Rational Mech. Anal. 130 
(1995), 343-400. 

[C2]. D. Christodoulou, Oral Communication (August 95). 

[CK]. D. Christodoulou and S. Klainerman, The Nonlinear Stability of the Minkowski space-time, Princeton Univ. 
Press, 1993. 

[CL]. D. Christodoulou and H. Lindblad, On the motion of the free surface of a liquid., Comm. Pure Appl. Math. 53 
(2000), 1536-1602. 

[CF]. R. Courant and K. O. Friedrichs, Supersonic flow and shock waves, Springer- Verlag, 1977. 

[Cr]. W. Craig, An existence theory for water waves and the Boussinesq and Korteweg-deVries scaling limits, Comm. 
in P. D. E. 10 (1985), 787-1003. 

[DM]. B. Dacorogna and J. Moser, On a partial differential equation involving the Jacobian determinant., Ann. Inst. 
H. Poincare Anal. Non. Lineaire 7 (1990), 1-26. 

[DN]. S. Dain and G. Nagy, Initial data for fluid bodies in general relativity, Phys. Rev. D 65 (2002). 

[El]. D. Ebin, The equations of motion of a perfect fluid with free boundary are not well posed., Comm. Part. Diff. 
Eq. 10 (1987), 1175-1201. 

[E2]. D. Ebin, Oral communication (November 1997). 

[F]. H. Friedrich, Evolution equations for gravitating ideal fluid bodies in general relativity, Phys. Rev. D 57 (1998). 

[FN]. H. Friedrich and G. Nagy, The initial boundary value problem for Einstein's vacuum field equation, Commmun. 
Math. Phys. 201 (1999), 619-655. 

[Ha]. R. Hamilton, Nash-Moser Inverse Function Theorem, Bull. Amer. Math. Soc. (N.S.) 7 (1982), 65-222. 

[HI]. Hormander, The boundary problem of Physical geodesy, Arch. Rational Mech. Anal. 62 (1976), 1—52.. 

[H2]. Hormander, Implicit function theorems (1977), Lecture Notes (Stanford ). 

[H3]. Hormander, The analysis of Linear Partial Differential Operators III, Springer Verlag. 

[Kl]. S. Klainerman, On the Nash-Moser-Hdrmander scheme, unpublished lecture notes. 

[K2]. S. Klainerman, Global solutions of nonlinar wave equations, Comm. Pure Appl. Math. 33 (1980), 43-101.. 

[LI]. H. Lindblad, Well posedness for the linearized motion of an incompressible liquid with free surface boundary., 
Comm. Pure Appl. Math., (2003). 

[L2]. , Well posedness for the linearized motion of a compressible liquid with free surface boundary., Comm. 

Math. Phys. (2003). 

[L3]. , Well posedness for the motion of an incompressible liquid with free surface boundary., preprint: 



http://xxx.lanl.gov/abs/math.AP/0402327, to appear in the Annals of Math.. 



[Na]. V.I. Nalimov,, The Cauchy-P oisson Problem (in Russian),, Dynamika Splosh. Sredy 18 (1974,), 104-210. 
[Ni]. T. Nishida, A note on a theorem of Nirenberg, J. Diff. Geometry 12 (1977), 629-633. 

[R]. A. D. Rendall, The initial value problem for a class of general relativistic fluid bodies, J. Math. Phys. (1992), 
1047-1053. 

[SY]. Schoen and Yau, Lectures on Differential Geometry, International Press, 1994. 

[S]. E. Stein, Singular Integrals and differentiability properties of functions (1970), Princeton University Press. 

[Wl]. S. Wu, Well-posedness in Sobolev spaces of the full water wave problem in 2-D, Invent. Math. 130 (1997), 39-72. 

[W2]. S. Wu, Well-posedness in Sobolev spaces of the full water wave problem in 3-D, J. Amer. Math. Soc. 12 (1999), 
445-495. 

[Y]. H. Yosihara, Gravity Waves on the Free Surface of an Incompressible Perfect Fluid 18 (1982), Publ. RIMS 
Kyoto Univ., 49-96. 



68 



